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Long-range corrected hybrid meta-generalized-gradient approximations with 

dispersion corrections 
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We propose a long-range corrected hybrid meta-GGA functional, based on a global 
hybrid meta-GGA functional, M05 [Y. Zhao, N. E. Schultz, and D. G. Truhlar, J. 



Chem. Phys. 123, 161103 (2005)], and empirical atom-atom dispersion correc- 
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tions. Our resulting functional, wM05-D, is shown to be accurate for a very wide 
range of applications, such as thermochemistry, kinetics, noncovalent interactions, 
equilibrium geometries, frontier orbital energies, fundamental gaps, and excitation 
energies. In addition, we present three new databases, IP131 (131 ionization poten- 
tials), EA115 (115 electron affinities), and FG115 (115 fundamental gaps), consisting 
of experimental molecular geometries and accurate reference values, which will be 
useful in the assessment of the accuracy of density functional approximations. 
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I. INTRODUCTION 



Because o 



(KS-DFT) 



[i 



its satisfactory accuracy and modest cost in many applications, Kohn-Sham density functional theory 
2j has become one of the most popular electronic structure methods for large ground-state systems 
3|-|6|. Its extension for treating excited-state systems, time-dependent density functional theory (TDDFT) {7, 8], 
has also been widely used. 

The crucial ingredient of KS-DFT, the exact exchange-correlation (XC) energy functional E xc [p], however, 
remains unknown and needs to be approximated. Functionals based on the local density approximation (LDA), 
modeling the XC energy densitylocally with that of a uniform electron gas (UEG), have been quite successful 
for nearly- free electron systems [3|, |4| , though still insufficiently accurate for most quantum chemical applications. 
Functionals based on the generalized gradient approximations (GGAs), additionally incorporating the gradient 
of local density into the LDA, have achieved reasonable accuracy in many applications. As an extension of the 
GGA (for rather restricted set of density variables), meta-GGA (MGGA) offers itself quite naturally. Functionals 
depending directly on the Laplacian of the density have not been pursued intensively, because of the difficulty of 
numerical evaluation. MGGAs, which adopt the kinetic energy density as a substitute for the Laplacian, have 



shown evidences of superiority over GGAs. (9 -111] 

However, the LDA, GGAs and MGGAs (commonly denoted as DFAs for density functional approximations) are 
based on the localized model XC holes, while the exact XC hole should be fully nonlocal. Currently, perhaps the 
most successful approaches to taking into account the nonlocality of XC hole are provided by hybrid DFT methods, 
incorporating a fraction of the exact Hartree-Fock (HF) exchange into the DFAs. Hybrid density functionals have 
achieved remarkable accuracy and have expanded the usefulness of DFT for many applications. Noticeably, global 



hybrid MGGA functionals 12l-l20|. where the XC energy density depends on the local density, the gradient of local 



density, a fraction of exact exchange, as well as the exact KS kinetic energy density (a function of the occupied KS 



orbitals) 



21 



24j . have been shown to potentially perform better than global hybrid GGA functionals 15H20l.l25ll26] 



due to the additional ingredient of kinetic energy density in global hybrid MGGA functionals. 

In global hybrid functionals, a small fraction of the exact HF exchange is added to a semilocal density functional. 
In certain situations, especially in the asymptotic regions of molecular systems, a large fraction (even 100%) of HF 
exchange is needed. Aiming to remedy this, long-range corrected (LC) hybrid DFT schemes have been actively 
developed 27H37j|. LC hybrids retain the full HF exchange only for the long-range electron-electron interactions, 
and thereby resolve a significant part of the self-interaction problems associated with global hybrid functionals. 

On the other hand, the development of accurate short-range (SR) exchange density functionals E x R [p], plays an 
important role in the progress of LC-DFT. In the first LC scheme, an ansatz for the conversion of any E x to E% R was 
proposed by Iikura et al. [271 ] . and has become widely used. However, their resulting LC hybrid GGA functionals 



3 



do not outperform the corresponding global hybrid GGA functional for thermochemistry. In 2006, Vydrov et al. 



proposed a different LC scheme 3JJ, based on integrating a GGA model exchange hole. Their resulting LC-wPBE 
functional has shown improved performance for thermochemistry and barrier heights, and is comparable to global 
hybrid GGA functionals such as B3LYP 



38 



. However, further improvements following this direction require 



the development of more accurate model exchange holes, which is a quite challenging task. 



Another approach to more accurate LC hybrid functionals was proposed by Chai and Head-Gordon 35|. First, 
augmenting the SR local spin density exchange energy density by a flexible enhancement factor (of the Becke's 



1997 form 



40j) and fully reoptimizing the LC functional on a diverse training set, yields the o;B97 functional. 



Second, including an adjustable fraction of SR HF exchange in the cuB97 functional with the similar reoptimization 
procedure, leads to the o;B97X functional. uB97 and wB97X have been shown to be accurate across a diverse set 
of test data, containing thermochemistry, kinetics, and noncovalent interactions [35[. 



However, problems associated with the lack of nonlocality of the DFA correlation hole, such as the lack of 
dispersion interactions (the missing of van der Waals forces), are not resolved by the LC hybrid schemes. The 
correlation functionals in typical LC hybrids are treated semilocally, which cannot capture the long-range (LR) 



correlation effects 



41 



42j . To remedy this, the DFT-D scheme was applied [43| to extend the cjB97X functional with 



damped atom-atom dispersion corrections, denoted as cjB97X-D [36|. Consequently, wB97X-D can obtain dispersive 



effects with essentially zero additional computational cost re 
has also been combined with the double-hybrid methods 44 



ative to cuB97X. As an alternative approach, wB97X 



481 ]. which mix both the HF exchange and nonlocal 



orbital correlation energy from the second-order perturbation energy expression in wave function theory. The 
resulting wB97X-2 functional [37] has yielded very high accuracy for thermochemistry, kinetics, and noncovalent 
interactions, though its fifth-order scaling with respect to system size may limit its applicability to larger systems. 



As the cjB97 series are LC hybrid GGAs, it seems a natural step to develop LC hybrid MGGAs and to assess 
their performance. In this work, we propose a new LC hybrid MGGA-D functional, denoted as wM05-D, which is 
shown to be accurate for a wide range of applications, when compared with the two closely related functionals: a 
global hybrid MGGA functional (M05-2X) [l9| and a LC hybrid GGA-D functional (wB97X-D) Q. The rest of 
this paper is organized as follows. In Sec. II, we briefly describe the relevant schemes developed in the LC hybrid 
approach. In Sec. Ill, we propose a new SR exchange functional, which serves as suitable basis functionals for 
systematically generating accurate LC hybrid MGGA functionals. The performance of the wM05-D functional is 
compared with other functionals in Sec. IV (on the training set), and in Sec. V (on some test sets). In Sec. VI, 
we give our conclusions. 
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II. RATIONALES OF LC HYBRID SCHEMES 

For the LC hybrid schemes, one first defines the long-range and short-range operators to partition the Coulomb 
operator. The most popular type of splitting operator used is the standard error function (erf), 

1 _ erf(o;ri2) + erfc(u;ri 2 ) ^ 

where r 12 = |r 12 | = |r x — r 2 | (atomic units are used throughout this paper). On the right hand side of Eq. ([I]), 
the first term is long-ranged, while the second term is short-ranged. The parameter u defines the range of these 
operators. 

In this work, we employ the erf/erfc partition, and use the following expression (as suggested in the recent LC 



hybrid schemes 



35 



36 



49|) for the LC hybrid functionals ( fractional number to be determined): 

£LC-DFA = ^LR-HF + ^SR-HF + (1 _ ^SR-DFA + £ DPA > (2) 

where E^ R ~ HF , the LR HF exchange, is computed by the occupied KS orbitals ^(r) with the LR operator, 

^ LR " HF = ~\ E E // ^(ri)^(r 2 )^^^(r 1 )^(r 2 )dr 1 dr 2 , (3) 

-Ef R ~ HF , the SR HF exchange, is computed similarly to the above but with the SR operator, 

^ SR " HF = -~ E E // C(ri)^(r 2 )^^^(r 1 )^(r 2 ) C /r 1 rfr 2 , (4) 

E^ R ~ DFA is the SR exchange approximated by DFAs , and E^ FA is the correlation functional the same as that of 
the full Coulomb interaction. 

In view of the E^~ BFA in Eq. 02]), as £"LR-hf anc j ^sr-hf are we rj (j e fi ne( j ; anc i accurate approximations for 
E^? FA are widely available, the accuracy of E^'^^ is thus closely related to the accuracy of a LC hybrid functional 



35 



361 ] . The analytical form of the SR LDA (the simplest SR DFA) exchange functional ^sr-lda can ^ e OD ^ a i nec i 



by the integration of the square of the LDA density matrix with the SR operator [5( 



Here, e^" LDA (p (T ) is the SR LDA exchange energy density for a-spin, 

e LDA (p.) = e^ A F(a a ), (6) 

where 

2 V4vr 



e^ A (p«) = 4(^) 1/3 pt /3 (r) (?) 
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is the LDA exchange energy density for cr-spin, kp a = (67r 2 p (T (r)) 1 ' /3 is the local Fermi wave vector, and a a 
u/{2kFa) is a dimensionless parameter controlling the value of the attenuation function F(a a ), 

n . . . . „ / i 



F(a a ) 



1 a c 

3 



A/vrerf 



2a„ 



3a CT + 4a£ + (2a CT - 4a£)exp 



4«L _ 



(8) 



To develop a possible SR DFA exchange functional E% R_DFA based on the knowledge of a DFA exchange functional 



E^ FA , there are three schemes as follows. Consider the general expression of DFA exchange functional, which is 



E 



DFA 



£ / e™ A (p a )F™ A dr, 



(9) 



F® FA is the DFA enhancement factor for u-spin. Depending on the type of DFA, F^J A = 1 for a LDA, F^J A 



Fg GA (p a , Vp ff ) for a GGA, F?J A = F^ GGA (p ff , Vp CT , r ff ) for a meta-GGA, where p CT (r) is the spin density, Vp CT (r) 
is the spin density gradient, and 



1 



7Vt 



(10) 



is the spin kinetic energy density. 

The first scheme was proposed by Iikura, Tsuneda, Yanai, and Hirao (ITYH) 
obtained by substituting a modified Fermi wave vector, 

kFa 



27 



28 



33|], where E^ R ~ DFA can be 



into SR exchange energy density of Eq. (Q, which a priori produces E^ R ~ DFA from any -E° FA , and reduces nicely 
to £?sr-lda from a ^lda Although the ITYH scheme ^possesses an admirable simplicity, some of its deficiencies 
(which potentially limit its accuracy) have been found 51|. 



The second scheme was proposed by Vydrov, Heyd, Krukau, and Scuseria (VHKS) [3l|, |32|, where for a given 
spherically-averaged exchange hole /i° FA (r, r 12 ), ^r-dfa - g evaluated as 



E 



SR-DFA 



2tt / p(r)dr / erfc^r^)/^ (r,r 12 )ri 2 dr 



12- 



(12) 



The pivot of this scheme is the engineering of the DFA exchange hole. The GGA model exchange hole of Ernzerhof 
and Perdew 52J (EP) provides a framework for modeling any GGA exchange hole. It has made considerable 
appearances in real applications after parametrization to reproduce the Perdew, Burke, and Ernzerhof (PBE) GGA 



[531 ]. In 2008, Henderson, Janesko, and Scuseria [5lJ (HJS) proposed another general model for the spherically 
averaged exchange hole corresponding to a GGA exchange functional, based on the work of EP. The HJS model 
improves upon the EP model by precisely reproducing the energy of the parent GGA, and by enabling fully analytic 
evaluation of range-separated hybrid density functionals. For meta-GGA, the TPSS exchange and correlation hole 



6 



models have been "reverse-engineered" 54|. However, the resulting LC-TPSS functional (a LC hybrid MGGA) has 



no satisfactory long-range correction effect 



3l| 



The third scheme was proposed by Chai and Head-Gordon (CHG) 35|, |36j, where E^ R ~ DFA is evaluated as 



E 



SR-DFA 



E 



SR-LDA 



(13) 



This simple scheme is expected to work well for a small u>. For highly parametrized E^ FA , such as the B97 [401 ]. 
M05 |l8j], and M08 26j functionals, the CHG scheme is particularly attractive due to its simplicity. However, 
how large is not too large for the u suitable for the CHG scheme? In the following sections, we will compare the 
performance of two new LC hybrid MGGA-D functionals, where one is developed by the CHG scheme, while the 
other is developed by a new scheme provided in this work, and our results help to answer the above question. 



III. LC HYBRID MGGA-D FUNCTIONALS 



In this section, we introduce our new LC hybrid MGGA-D functionals. Note that LC-TPSS has been developed 
by utilizing the TPSS exchange hole (based on the VHKS scheme) 



31], but it is found that LC-TPSS does not 
benefit much by admixture of HF exchange. The Mil functional 49] has been developed based on the extension 



of a global hybrid MGGA functional, M08 
Parallel to the strategy of the cjB97 series 



26|Lto 



35 



35] 



C-DFT, following the CHG scheme 

, we choose to modify the M05 functional. The M05 functional 
is a global hybrid MGGA functional with a powerful form 18j, Il9| . and our work is based on modifying this 
functional. Its exchange part consists of the PBE exchange functional and a reasonable kinetic-energy-density 
enhancement factor. The PBE exchange is a theoretically sound starting point because it satisfies the correct UEG 
limit and also has reasonable behavior at large values of the reduced spin density gradient s a . 

To satisfy the UEG limit of SR exchange, we replace the PBE exchange energy density e PBE (p CT , Vp CT ) with 
the SR-PBE exchange energy density e^" PBE (p (J , Vp CT ) generated by the HJS model exchange hole (based on the 
VHKS scheme), whose virtues are indicated in Sec. II. To achieve a flexible functional form, we retain the kinetic- 
energy- density enhancement factor (similar to the CHG scheme). We denote this resulting functional as SR-M05 
(short-range M05) exchange, as it reduces to the M05 exchange at u — 0. 



(14) 



£SR-M05 = £ / e PBE ( P ., Vp a )f(w a )dT, 

cr J 

where f{w a ) is the kinetic-energy-density enhancement factor, 

in 

i=0 



(15) 
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w a is a function of i CT , and t a is a function of the kinetic energy density r CT of electrons with spin a, as designed by 



Becke 



where 



23|, 



w a = (U - l)/(t a + 1), 



r LDA ^ ^ (67r2)2 /3 p 5/3_ 



(16) 

(17) 
(18) 



In general, the enhancement factor should be w-dependent. But from the works of LC-TPSS 3lJ and Mil |49j], the 
optimal bj for a LC hybrid MGGA is expected to be small as well. For a sufficiently small u value, our proposed 
functional form, inspired by the VHKS and CHG schemes, should be a good approximation. 

We use the same form for the correlation functional as the M05 correlation functional, which can be decomposed 
into same-spin E^ 5 a and opposite-spin E^ 1 ^ components, 



^M05 _ ^M05 _|_ 



M05 
c a/3' 



For the opposite-spin terms, 



F M05 _ / LDA / 2 \ t 
a/3 — / e c a^9a^ av )ar, 



fa/3 («o 



i=Q 



U 



a/3 



1 + 7a/JS^ 



and for the same-spin terms, 



7q/3 = 0.0062, 

2 1 / 2 , 2 \ 

Sail _ 2 S ' 9 

/W 
^ w (4)(l-^-)dr, 
To- 

n 

i=0 
lcra S cr 



1 + la* S 2 a 
lea = 0.06. 



(19) 

(20) 
(21) 

(22) 
(23) 
(24) 

(25) 
(26) 

(27) 
(28) 



1 — tJ^/to- is a sel 
energy density 



interaction correction factor proposed by Becke, [22] in which is the von Weizsacker kinetic 
55| given by 

w l v P<x| 2 



(29) 
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In a one-electron case, r a = , so Eq. (|25|) vanishes in any one-electron system. The correlation e nerg y densities 
e c I Q^ an d e c I o-a : are derived from the Perdew-Wang parametrization of the LDA correlation energy 56|, using the 
approach of Stoll et al. 57| . 



3 LDA 

'c a/3 



,LDA/ 



(p a ,P^)-er A (p a ,0)-e^ A (0,p,) 



e^ = e^ DA (p CT ,0). 



(30) 



(31) 



Based on the above functional expansions, we propose a new LC hybrid MGGA functional, wM05-D. It contains 
a fraction of the SR HF exchange, 



xjM05-D 



^LR-HF _|_ c ^SK-Ht _|_ ^s±t-MUt> _j_ 

02 32 £Q 



SR-HF , t-iSR-MOS 



nM05 



(32) 



We enforce the exact UEG limit for the wM05-D functional by imposing the following constraints: 

C CT<X,0 = 1) 



and 



c a/3,o — 1> 



a + c x = 1. 



Following the general form of the DFT-D scheme [43j, our total energy 



DFT-D — -^KS-DFT "1" -C/disp 



+ E d 



(33) 



(34) 



(35) 



(36) 



is computed as the sum of a KS-DFT part and an empirical atomic-pairwise dispersion correction. We choose to 
use the same form of unsealed dispersion correction as implemented in co>B97X-D 361 ]. 

ATat-l 



-Edisp — ~~ ^ d ^ d "^g"/damp(-Rjj) ? 
i=l j=i+l U ij 



(37) 



where N at is the number of atoms in the system, C$ is the dispersion coefficient for atom pair ij, and Rij is an 
interatomic distance. The damping function, 

1 



/damp(^)- l + a ( R ../ Rr) -12 



(38) 



enforces the conditions of zero dispersion correction at short interatomic separations and correct asymptotic pairwise 
vdW potentials. Here, R r is the sum of vdW radii of the atomic pair ij, and the only non-linear parameter, a, 
controls the strength of dispersion corrections. 
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To achieve an optimized functional 



59 



or well-balanced performance across typical applications, we use the same 
351 ]. which contains 412 accurate experimental and accurate theoretical results, 

], the atomization energies of the G3/99 set 
61] . the ionization potentials (IPs) of the G2-1 set 62] (40 molecules, excluding SH 2 ( 2 A 1 ) and 



diverse training set described in Ref. 

including the 18 atomic energies from the H atom to the Ar atom 
(223 molecules 



N 2 ( 2 I1) cations due to the known convergence problems for semilocal density functionals |60j), the electron affinities 
(EAs) of the G2-1 set (25 molecules), the proton affinities (PAs) of the G2-1 set (8 molecules), the 76 barrier heights 



15 



of the NHTBH38/04 and HTBH38/04 sets 
data is weighted ten times more than the others. A 
by a least-square fitting procedure described in Ref. 



63], and the 22 noncovalent interactions of the S22 set 64|. The S22 
1 the parameters in co>M05-D are determined self-consistently 
35]. For the non-linear parameter optimization, we focus on a 



range of possible u values (0.0, 0.1, 0.2, 0.3, and 0.4 Bohr x ), and optimize the corresponding a values in the steps 



described in Ref. 



36]. 



M05 and M05-2X 



18 



191 ] both used m=ll in Eq. f|15p and n=4 in Eqs. ( 12 ip and ( I26p . However, during the 



optimization procedure of wM05-D, we found that the statistical errors are close for m=10 and m=ll, while the 
one with m=ll has parameters significantly larger. A recent study by Wheeler and Houk has shown that large 
magnitude of the parameters in Eq. ( TT5]) may result in large grid errors 65] . Moreover, the use of large parameters 
increases the possibility of convergence difficulty as well as the over-fitting effects. Thus, we choose m=10 instead 
of 11 in Eq. (175]) . The optimized parameters of the o;M05-D functional are given in Table [IJ in which the u> value 
is same as that of wB97X-D, while the fraction of SR HF exchange, c x , is larger than that of wB97X-D (^ 0.22). 
This helps to reduce the self-interaction error (SIE) of the functional, as can be seen in Sec. V. 

We also tried a simple model (based on the CHG scheme), where the SR-PBE exchange energy density 
e^" PBE (p CT , Vp CT ) used in Eq. (IH]) is substituted with e^" LDA (p (T )F ;c PBE (s ff ), that is, the SR LDA exchange en- 
ergy density in Eq. (JH]) multiplied by the PBE enhancement factor. We tried this because the mathematical form 
of the latter is significantly simpler than that of the former, and is the model on which Mil based. The parametriza- 
tion is the same for this simple model, which we denoted by o;M05s-D. Compared to wM05-D, the optimal u value 
is also 0.2 bohr -1 , but the corresponding optimal a value is found to be 100 and the linear parameters are also 
larger. 



IV. RESULTS FOR THE TRAINING SET 



All calculations are performed with a development version of Q-CHEM 3.2 [66]. Spin-restricted theory is used 
for singlet states and spin-unrestricted theory for others, unless noted otherwise. For the binding energies of the 



weakly bound systems, the counterpoise correction 



671 ] is employed to reduce basis set superposition error (BSSE). 



Results for the training set are computed using the 6-311++G(3df,3pd) basis set with the fine grid, EML(75,302), 
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TABLE I. Optimized parameters for wM05-D. Here, the non-linear parameter a is denned in Eq. (J38]), and others are 
defined in Eq. fl32]) 



a 30.0 

lo 0.2 Bohr" 1 

c x 0.369592 

0.630408 1.00000 1.00000 

1 -0.219121 -0.95491 -5.26863 

2 -0.14411 12.138 17.9935 

3 1.27732 -35.1041 -17.6408 

4 -1.59959 19.5804 0.625687 

5 -5.94702 

6 13.5822 

7 10.5048 

8 -28.7168 

9 -6.89761 

10 19.0574 



consisting of 75 Euler-Maclaurin radial grid points [68] and 302 Lebedev angular grid points [69]. The error for 
each entry is defined as error = theoretical value — reference value. The notation used for characterizing statistical 
errors is as follows: mean signed errors (MSEs), mean absolute errors (MAEs), root-mean-square (rms) errors, 
maximum negative errors (Max(— )), and maximum positive errors (Max(+)). 

First, we show the results of the first iteration of fitting procedure, comparing the new LC scheme with the 
CHG scheme (the simple model) for o;=0.1, 0.2, 0.3 and 0.4 bohr -1 . We optimize o;M05 and wM05s using the 
corresponding wPBE and cuPBEs orbitals, and denote these optimized functionals as wM05* and cjM05s*. The 
statistical errors are believed to be quite close to those obtained self-consistently. As can be seen in Table [TlJ the 
difference between the performance of cuM05* and wM05s* is noticeable for w=0.2 bohr -1 , and becomes larger for a 
larger u value. Therefore, a LC hybrid MGGA functional with a larger u value (such as Mil with w=0.25 bohr -1 ) 
may perform better with our new scheme than with the CHG scheme. 

In subsequent iterations, we include the dispersion corrections, increase the training weight of S22 set, and found 
the functionals optimized with u=0.2 bohr -1 . To view the effect of the long-range correction and the dispersion 



11 



corrections, we also consider the functional form M05 and M05-D. The latter is the limiting case where oo=0 for 
cjM05-D, of which the corresponding optimal a value is found to be 2. We reoptimize M05 and M05-D functionals 
on the same training set using the M05-2X orbitals, truncate their functional expansions at the same orders m=10 
and n=4, and denote these two reoptimized functionals as M05* and M05-D*. Just like the o;B97X functional 
without dispersion correction, all data in the training set are equally weighted in the least-squares fitting for M05*. 



The overal 



performance of our new cjM05-D is compared with the trial simple model cuM05s-D, M05-D*, M05* 



and M05-2X [19f, as well as existing wB97X-D (a LC hybrid GGA-D) [36|. Note that M05 [18] and M05-2X share 
the same functional form, but the former is distracted to deal with transition-metal compounds, so the latter should 
be our concern. In the o;B97 series, cuB97X-D has the closest relationship to u;M05-D, while o;B97 and cjB97X, 
developed without dispersion corrections, are expected to perform poorly for noncovalent interactions. 

In Table IHIt the first comparison (a;M05-D vs. wM05s-D) partially determines the choice of our proposed 
functional. Although u;M05-D performs worse than co>M05s-D for HTBH, the overall performance of cjM05-D in 
the training set is the best. 

A second comparison between co>M05-D and M05-D* indicates that the exact long-r ang e exchange indeed leads 



31 



35 



The third comparison 



to an overall improvement to MGGA, although not as large as that to GGA 
is between M05-D* and M05*. The cooperation of the training weight and the empirical dispersion corrections 
leads to a significant improvement in the results for noncovalent interactions (the S22 data) and a modest overall 



change. Recently, there have been the updated reference values for the S22 set 70]. We have also examined the 



7l| 



performance of wM05-D against the updated S22 reference values. As shown in the supplementary material 
the overall performance of the functional against the updated reference values is similar to that against the original 
ones. 



V. RESULTS FOR THE TEST SETS 



To test the performance of cjM05-D outside its training set, we also evaluate its performance on various test 
sets involving 48 atomization energies in the G3/05 test set 72|, 30 chemical reaction energies taken from the 



15 



63(, 29 noncovalent interactions 



63 



64|, 166 optimized geometry 



NHTBH38/04 and HTBH38/04 databases 
properties of covalent systems [731 ]. 12 intermolecular bond lengths 64|, 4 dissociation curves of symmetric radical 



cations as well as three new databases, consisting of 131 vertical IPs, 115 vertical EAs and 115 fundamental gaps. 
For excitation energies, we perform TDDFT calculations for 19 valence excitation energies, 23 Rydberg excitation 
energies and one long-range charge transfer excitation curve of two well-separated molecules. Each EA can be 
evaluated by two different ways, and each fundamental gap can be evaluated by three different ways, so there are 
a total of 1038 pieces of data in the test sets, which are larger and more diverse than the training set. Unspecified 
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TABLE II. Comparisons between the u;M05* and a;M05s* functionals (denned in the text) for different co values. Statistical 
errors are in kcal/mol. 



System 


u (Bohr" 1 ) 
Error 


0.1 

uM05* 


0.1 

wM05s* 


0.2 

wM05* 


0.2 

wM05s* 


0.3 

wM05* 


0.3 

wM05s* 


0.4 

uM05* 


0.4 

wM05s* 


Atoms 


MSE 


-0.15 


0.24 


0.05 


0.63 


0.23 


0.97 


0.46 


1.33 


(18) 


MAE 


2.02 


2.09 


1.81 


2.35 


2.00 


3.36 


3.22 


5.05 


G3/99 


MSE 


0.06 


0.09 


0.05 


0.03 


-0.04 


-0.12 


-0.18 


-0.27 


(223) 


MAE 


1.77 


1.79 


1.66 


1.76 


1.78 


2.02 


2.10 


2.35 


IP 


MSE 


-0.58 


-1.48 


-0.84 


-1.36 


-0.38 


-0.32 


0.30 


0.73 


(40) 


MAE 


2.75 


3.06 


2.81 


3.08 


2.68 


2.81 


2.64 


2.79 


EA 


MSE 


-1.50 


-1.70 


-1.29 


-1.15 


-0.94 


-0.70 


-0.64 


-0.39 


(25) 


MAE 


2.50 


2.56 


2.33 


2.22 


2.07 


1.97 


1.98 


1.91 


PA 


MSE 


-1.65 


-2.68 


-1.49 


-2.71 


-1.07 


-2.11 


-0.78 


-1.54 


(8) 


MAE 


1.87 


2.68 


1.83 


2.82 


1.79 


2.53 


1.86 


2.43 


NHTBH 


MSE 


-1.26 


-1.09 


-0.68 


-0.39 


0.08 


0.40 


0.85 


1.17 


(38) 


MAE 


1.98 


1.82 


1.51 


1.40 


1.46 


1.67 


1.71 


1.95 


HTBH 


MSE 


-1.96 


-1.95 


-1.95 


-1.68 


-1.61 


-1.24 


-1.21 


-0.84 


(38) 


MAE 


2.19 


2.08 


2.12 


1.84 


1.95 


1.56 


1.86 


1.47 


S22 


MSE 


2.65 


1.91 


1.80 


1.06 


1.01 


0.56 


0.45 


0.26 


(22) 


MAE 


2.65 


1.91 


1.80 


1.06 


1.02 


0.67 


0.71 


0.63 


All 


MSE 


-0.31 


-0.42 


-0.30 


-0.35 


-0.21 


-0.18 


-0.11 


-0.02 


(412) 


MAE 


2.03 


2.03 


1.86 


1.90 


1.84 


2.01 


2.06 


2.28 
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TABLE III. Statistical errors (in kcal/mol) of the training set. The M05-D* and M05* functionals are defined in the text. 
M05-2X was not particularly parametrized using this training set. 



System 


Error 


wM05-D 


wM05s-D 


M05-D* 


M05* 


M05-2X 


wB97X-D 


Atoms 


MSE 


0.37 


0.83 


0.18 


-0.48 


-3.01 


-0.05 


(18) 


MAE 


2.02 


2.28 


2.61 


1.98 


5.10 


2.57 


G3/99 


MSE 


-0.03 


-0.03 


-0.10 


-0.05 


2.01 


-0.24 


(223) 


MAE 


1.62 


1.73 


1.78 


1.78 


3.65 


1.93 


IP 


MSE 


-0.80 


-1.33 


0.06 


0.27 


1.10 


0.19 


(40) 


MAE 


2.86 


3.04 


2.84 


2.51 


3.35 


2.74 


EA 


MSE 


-1.02 


-0.98 


-0.54 


-0.84 


-0.23 


0.07 


(25) 


MAE 


2.12 


2.13 


2.13 


2.35 


2.48 


1.91 


PA 


MSE 


-1.48 


-2.66 


-0.94 


-1.76 


-1.26 


1.42 


(8) 


MAE 


2.10 


3.07 


1.31 


2.17 


1.51 


1.50 


NHTBH 


MSE 


-0.94 


-0.59 


-1.38 


-1.32 


0.13 


-0.45 


(38) 


MAE 


1.57 


1.53 


2.04 


2.08 


1.75 


1.51 


HTBH 


MSE 


-2.82 


-2.33 


-2.95 


-1.77 


-0.65 


-2.57 


(38) 


MAE 


2.83 


2.37 


3.08 


2.14 


1.51 


2.70 


S22 


MSE 


-0.01 


-0.01 


0.04 


3.46 


0.73 


-0.08 


(22) 


MAE 


0.27 


0.21 


0.23 


3.46 


0.87 


0.21 


All 


MSE 


-0.51 


-0.49 


-0.49 


-0.21 


1.02 


-0.36 


(412) 


MAE 


1.83 


1.89 


1.99 


2.05 


3.05 


1.96 



14 



TABLE IV. Statistical errors (in kcal/mol) of the test sets. 



System Error wM05-D wM05s-D M05-2X wB97X-D 

G3/05 MSE -0.85 -1.67 0.00 0.25 

(48) MAE 3.21 3.79 5.24 3.02 

RE MSE -0.58 -0.65 -0.86 -0.24 

(30) MAE 1.49 1.32 1.65 1.63 

Non-covalent MSE -0.11 -0.05 0.50 -0.15 

(29) MAE 0.31 0.30 0.61 0.43 

All MSE -0.58 -0.95 -0.11 0.01 

(107) MAE 1.94 2.15 2.98 1.93 

detailed information of the test sets as well as the basis sets, and numerical grids used is given in Ref . |35| . 

A. Atomization Energies, Reaction Energies, and Noncovalent Interactions 



Table IIVI summarized the general energetic results in the same way as in Ref. 36] , for convenience of further 
comparisons. Since the 30 chemical reaction energies are taken from the NHTBH38/04 and HTBH38/04 databases 
calculated in Table IITT| the EML(75,302) grid is used. In Table HVt the comparison between cjM05-D and wM05s- 
D shows noticeable difference in atomization energies, and makes great influence on the choice of our proposed 
functional. 



B. Equilibrium Geometries 



Satisfactory predictions of molecular geometries of covalent and non-covalent systems by density functionals 
are necessary for practical use. For covalent systems, we perform geometry optimizations for each functional on 

n 

the equilibrium experimental test set (EXTS) [73(, while for non-covalent systems, we compute the intermolecular 
bond lengths of 12 weakly bound complexes taken from the S22 set [64], using 6-311+- |-G(3df,3pd) basis set with 
the EML(75,302) grid. As shown in Table |V] performance of all the hybrid functionals in predicting optimized 
geometries of EXTS is similar, while the performance of simple model (a;M05s-D) is somewhat worse for the 
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TABLE V. Statistical errors (in A) of EXTS (Ref. [73J]) and bond lengths of 12 weakly bound complexes from the S22 set 



(Ref. [64]). The results of wB97X-D are taken from Ref. 



36]. 



System 



Error 



wM05-D 



wM05s-D 



M05-2X 



wB97X-D 





MSE 


0.003 


0.001 


-0.004 


-0.002 




MAE 


0.010 


0.009 


0.009 


0.009 


EXTS (166) 


rms 


0.019 


0.014 


0.014 


0.013 




Max(-) 


-0.081 


-0.083 


-0.082 


-0.078 




Max(+) 


0.177 


0.067 


0.054 


0.055 




MSE 


-0.041 


-0.069 


-0.021 


-0.044 




MAE 


0.061 


0.078 


0.062 


0.064 


Weak (12) 


rms 


0.083 


0.102 


0.080 


0.085 




Max(-) 


-0.189 


-0.195 


-0.165 


-0.198 




Max(+) 


0.043 


0.029 


0.140 


0.056 



intermolecular bond lengths. We decide our proposed model to be wM05-D in this subsection. For brevity, the 
performance of o;M05s-D will not be shown for subsequent calculations. 



C. Dissociation of Symmetric Radical Cations 



Common semilocal functionals are generally accurate for systems near equilibrium. However, due to considerable 



self-interaction errors in semilocal functionals, spurious fractional charge dissociation occurs 



mm- This 



situation becomes amplified for symmetric charged radicals Xj, such as Hj", He^~, Ne^~ and ArJ. Grafenstein and 



coworkers have obtained qualitatively correct result for these systems 



76 



771 ] using self-interaction-corrected DFT 



proposed by Perdew and Zunger 78|, and confirmed that the errors of standard DFT methods should be dominated 
by the SIEs. 

We perform unrestricted calculations with the aug-cc-pVQZ basis set and a high-quality EML (250,590) grid. 
The DFT results are compared with results from HF theory, and the very accurate CCSD(T) theory 79j, [80J. The 
HF method is exact in Fig. (TJ and gives qualitatively correct results from Fig. [2] to Fig. HI Although wM05-D has 
the same amount of LR HF exchange as wB97X-D, the larger fraction of SR HF exchange included in wM05-D 
helps to reduce its remaining SIE. Therefore, the error of wM05-D is smaller than that of wB97X-D, especially for 
larger cations (e.g. Ne^~ and Ar^~). The global hybrid functional M05-2X exhibits the undesirable X^~ dissociation 
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FIG. 1. Dissociation curve of H^. Zero level is set to E(H)+E(H + ) for each method, 
curves, displaying a spurious energy barrier at intermediate bond length R. 



D. Frontier Orbital Energies 

Let IP (AT) be the ionization potential and EA(iV) be the electron affinity of the A^-electron system, which are 
defined as 

IP(N) = E N ^-E N , (39) 

EA(iV) = E N - E N+1 , (40) 

respectively, with being the total energy of A^-electron system. For the exact DFT, the vertical ionization 
potential of a neutral molecule is identical to the minus HOMO (highest occupied molecular orbital) energy of the 



neutral molecule 



fill 



IP(iV) = -e N (N), (41) 

and the vertical electron affinity of a neutral molecule is identical to the minus HOMO energy of the anion (since 
EA(A0 = IP (AT + 1) by definition), 

EA(N) = -e N+1 (N + l), (42) 
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FIG. 2. Dissociation curve of He^. Zero level is set to E(He)+E(He + ) for each method. 



where 6m (N) is the M-th orbital energy of iV-electron system. The vertical electron affinity of a neutral molecule 
may also be approximated by the minus LUMO (lowest unoccupied molecular orbital) energy of the neutral 
molecule, but it is proved that there exists a difference between the vertical EA and the minus LUMO energy, 



A xc = e N+1 (N+l)-e N+1 {N), 



(43) 



where the difference A xc arises from the discontinuity of exchange-correlation potentials 8214841] . Recent study 



85], so the minus LUMO energy calculated by a LC 



shows that A xc is close to zero for LC hybrid functionals 
hybrid functional should be close to the vertical EA. 

To evaluate the performance of the functionals on the HOMO energy of the neutral molecule, we collect a new 
database, IP131, which consists of experimental vertical IPs of 18 atoms and 113 molecules in the experimental 



geometries, 
publications 



he geometries and most of the reference values are collected from the NIST database [86j. Other 



871 ] are adopted for the experimental vertical IPs of some molecules. The DFT calculations are 



performed with 6-311++G(3df,3pd) basis and EML(75,302) grid. As can be seen in Table [VTl o;M05-D gives the 
best results. The global hybrid M05-2X gives the worst results here, due to its incorrect long-range XC-potential 
behavior. 

To evaluate the performance of the functionals on the vertical electron affinity, we construct another database 
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FIG. 3. Dissociation curve of Nejj - . Zero level is set to E(Ne)+E(Ne + ) for each method. 



called EA115, which consists of 18 atoms and 97 molecules. For the molecular geometries, it is a subset of IP131. 
Because experimental vertical EAs are not as widely available as experimental vertical IPs, the reference values 
of vertical EAs are obtained via the accurate CCSD(T) calculations (using Eq. (|40|) ). The CCSD(T) correlation 
ener gies in the basis-set limit are extrapolated from calculations using the aug-cc-pVTZ and aug-cc-pVQZ basis 
sets 



picorr v3 rpcoTT\^3 
- x ~ Y (AA\ 



where X=3 and Y =4 for the aug-cc-pVTZ and aug-cc-pVQZ basis, respectively. The electron affinities are evaluated 
in two different ways, as shown in Table IVDl for the minus HOMO energy of the anion, and Table IVTHi for the minus 
LUMO energy of the neutral molecule. Clearly, the LC hybrid functionals outperform the global hybrid M05-2X. 



The reference values and molecular geometries of IP131 and EA115 are given in the supplementary material 
along with detailed DFT results. 
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-50 



10 



R (angstrom) 



FIG. 4. Dissociation curve of Ar^. Zero level is set to E(Ar)+E(Ar + ) for each method. 



TABLE VI. Statistical errors (in eV) for the IP131 database. Error is defined as - €n(N) - IP ver ti ca i • Experimental geometries 
and reference values are used for all molecules. 

System Error wM05-D M05-2X wB97X-D 

atoms MSE -1.48 -2.06 -1.64 

(18) MAE 1.48 2.06 1.64 

rms 1.74 2.16 1.98 

molecules MSE -0.68 -1.23 -0.92 

(113) MAE 0.68 1.23 0.92 

rms 0.76 1.27 1.00 

total MSE -0.79 -1.34 -1.02 

(131) MAE 0.79 1.34 1.02 

rms 0.96 1.43 1.18 
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TABLE VII. Statistical errors (in eV) for the EA115 database. Error is defined as - €n+i(N + 1) - EA vcrt i ca i. Experimental 
geometries and CCSD(T) reference values are used for all molecules. 

System Error wM05-D M05-2X wB97X-D 

atoms MSE -0.46 -1.21 -0.53 

(18) MAE 0.49 1.21 0.57 

rms 0.73 1.35 0.84 

moelcules MSE -0.54 -1.18 -0.54 

(97) MAE 0.55 1.18 0.56 

rms 0.80 1.32 0.82 

total MSE -0.53 -1.18 -0.54 

(115) MAE 0.55 1.18 0.56 

rms 0.79 1.32 0.82 



TABLE VIII. Statistical errors (in eV) of the minus LUMO energy of the neutral molecule for the EA115 database. 



Experimental geometries and CCSD(T) reference values are used for all molecules. 

System Error wM05-D M05-2X wB97X-D 

atoms MSE -0.27 0.57 -0.02 

(18) MAE 0.73 1.02 0.74 

rms 0.92 1.12 0.89 

moelcules MSE -0.24 0.60 0.05 

(97) MAE 0.60 0.75 0.52 

rms 0.69 0.94 0.60 

total MSE -0.24 0.60 0.04 

(115) MAE 0.62 0.79 0.55 

rms 0.73 0.97 0.65 



E. Fundamental Gaps 

The fundamental gap E g of a molecule with N electrons is defined as 

E g = IP(N) -EA(JV). (45) 
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Following Eqs. ( 139|) and (1401) for the definitions of IP and EA, three self-consistent field (SCF) calculations (for the 
neutral molecule, cation and anion) are required to obtain the fundamental gap of a molecule. Using Eq. (14 ip and 
( |4"2"1) . the fundamental gap of a molecule can also be obtained by two SCF calculations (for the neutral molecule 
and anion). 

Following Janak's theorem 



the fundamental gap can be approximated by the HOMO-LUMO gap 



84] 



A KS = e N+1 (N)-e N (N), (46) 



and we can obtain the fundamental gap of a system using only one calculation. But from Eqs. ( 14TI) . ( l4~2l) . ( )43l) . 
f|45l) . and f|46l) . we know that there exists a difference between the fundamental gap and HOMO-LUMO gap, 

E g = A KS + A xc . (47) 

so the HOMO- 



As previously mentioned, A xc has been shown to be close to zero for LC hybrid functionals 
LUMO gap calculated by a LC hybrid functional should be close to the fundamental gap. 

To evaluate the performance of the functionals on fundamental gap, we construct another database called 
FG115, which shares the same molecular geometries with the EA115 database. For consistency, the reference 
values of fundamental gaps are also obtained via the CCSD(T) calculations described in the last subsection (using 
Eqs. (1391), gO]), and (pE}). 

To examine the performance of density functionals, we evaluate the fundamental gaps using three different 
estimates, with 6-311++G(3df,3pd) basis and EML(75,302) grid. The results are shown from Table HXl to Table 
IXH in order of increasing the number of SCF calculations required for each molecule. In the estimate requiring three 
calculations, the results are similar for the three functionals. wB97X-D gives worse results than other functionals 
in the estimate requiring two calculations. In the simplest estimate, the HOMO-LUMO gap, which requires only 
one SCF calculation for each system, cuM05-D significantly outperforms the other two functionals. The reference 
values of FG115 and detailed HOMO-LUMO gap results by DFT methods are given in the supplementary material 
HI 



F. Excitation Energies 



To assess the performance of density functionals on excitation energies, we perform TDDFT calculations on five 



small molecules 



which include nitrogen gas (N 2 ), carbon monoxide (CO), water (H 2 0), ethylene (C 2 H 4 ) and 



formaldehyde (CH 2 0), with 6-311(2+, 2+)G** basis and EML(99,590) grid. The molecular geometries, experimen- 



tal values of excitation energy are taken from Ref. 90j. The detail results and mean absolute errors of all excited 
states are listed in Table Kill The new cjM05-D functional yields excellent performance, especially for the Rydberg 
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TABLE IX. Statistic errors (in eV) of HOMO-LUMO gaps for the FG115 database. The energy gap of each system is 



evaluated by only one SCF calculation. 

System Error wM05-D M05-2X u;B97X-D 

atoms MSE -1.14 -2.56 -1.55 

(18) MAE 1.43 2.56 1.79 

rms 1.62 2.79 2.05 

molecules MSE -0.62 -2.00 -1.15 

(97) MAE 0.73 2.00 1.15 

rms 0.93 2.13 1.34 

total MSE -0.70 -2.08 -1.21 

(115) MAE 0.84 2.08 1.25 

rms 1.07 2.24 1.48 



TABLE X. Statistic errors (in eV) of fundamental gaps for the FG115 database, each evaluated by the difference of HOMO 
energies between the neutral molecule and anion. The energy gap of each system is evaluated by two SCF calculations. 



System Error wM05-D M05-2X wB97X-D 

atoms MSE -0.95 -0.83 -1.04 

(18) MAE 0.98 0.87 1.08 

rms 1.17 1.00 1.30 

molecules MSE -0.31 -0.42 -0.55 

(97) MAE 0.56 0.51 0.72 

rms 0.70 0.60 0.85 

total MSE -0.41 -0.48 -0.63 

(115) MAE 0.62 0.57 0.78 

rms 0.79 0.68 0.93 



excitations. Note that wM05-D outperforms wB97X-D in both HOMO energies and Rydberg excitations, due to 
the larger fraction of short-range HF exchange included in wM05-D (both functionals possess the same amount of 
LR-HF exchange). 

Following Dreuw et al, we perform TDDFT calculations for the lowest charge-transfer (CT) excitation between 
ethylene and tetrafluroethylene, with a separation of R. Dreuw et al. have shown that the exact CT excitation 
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TABLE XL Statistic errors (in eV) of IP-EA values for the FG115 database. The energy gap of each system is evaluated 



by three SCF calculations. 

System Error wM05-D M05-2X o;B97X-D 

atoms MSE 0.28 0.33 0.28 

(18) MAE 0.35 0.36 0.36 

rms 0.60 0.63 0.59 

molecules MSE 0.34 0.42 0.22 

(97) MAE 0.44 0.50 0.39 

rms 0.73 0.78 0.68 

total MSE 0.33 0.40 0.23 

(115) MAE 0.43 0.48 0.39 

rms 0.71 0.75 0.66 



energy from the HOMO of donor to the LUMO of acceptor should have the following asymptote 9l| : 



R 



IPt 



EA/ 



(48) 



where IPd is the ionization potential of donor and EAa is the electron affinity of acceptor. Fig. |5] shows the trend of 
the excitation curves, and indicates the LC hybrid functionals obviously outperforms the global hybrid M05-2X. For 
the values of the excitation energies, as shown in Fig. El u;M05-D is about 0.2 electron volt better than wB97X-D. 



VI. CONCLUSIONS 



We have developed a LC hybrid MGGA-D functional, called wM05-D, which includes 100% long-range exact 
exchange, a fraction (« 37%) of short-range exact exchange, a modified M05 exchange density functional for 



short-range interaction, the M05 correlation density functional [18|, [19|], and empirical atomic-pairwise dispersion 
corrections. For the modified short-range M05 exchange density functional, we have investigated two models. After 
comparisons in the training set and test sets, we decide to propose the one based on our new LC scheme cuM05-D, 
and marked the trial one as co>M05s-D. When the constraint of u = is applied, a;M05-D and wM05s-D are both 



reduced to the existing M05 functional form 



18. 



19] with the same empirical atomic-pairwise dispersion corrections. 



The constrained form (u = 0), when re-optimized on the same training set, provides worse performance on the 
training set, indicating that the single extra degree of freedom corresponding to long-range exchange is of physical 
significance to a hybrid MGGA. 
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TABLE XII. Vertical excitation energies (in eV) of several low- lying excited states of N2, CO, water, formaldehyde and 



ethylene using 6-311(2+, 2+)G** basis set. The geometries and experimental values are taken from Ref. [90' 



Mol. 


State 


Exp. 


wM05-D 


M05-2X 


wB97X-D 




v 1 n 9 


9.31 


9.30 


9.42 


9.38 






9.97 


8.76 


8.35 


9.31 






10.27 


10.14 


10.51 


9.82 


N 2 


V 3 £+ 


7.75 


7.86 


8.30 


7.17 




v 3 n 9 


8.04 


7.94 


8.12 


7.82 




V 3 A U 


8.88 


8.74 


8.35 


8.23 




V 3 E" 


9.67 


8.76 


9.26 


9.31 




v 3 n„ 


11.19 


11.30 


11.72 


10.98 




v x n 


8.51 


8.51 


8.74 


8.47 






9.88 


9.36 


9.11 


9.78 


CO 


v 3 n 


6.32 


6.66 


7.03 


6.07 




v 3 s+ 


8.51 


8.47 


8.87 


8.00 




V 3 A 


9.36 


9.19 


9.11 


8.88 




V 3 E" 


9.88 


9.36 


9.55 


9.78 




R x Bi 


7.4 


7.68 


8.04 


7.23 




R 1 A 2 


9.1 


9.14 


9.60 


8.63 


H 2 


R 1 Ai 


9.7 


9.73 


10.29 


9.20 




R J Bi 


10.0 


9.72 


10.32 


9.17 




R x Ai 


10.17 


10.06 


10.71 


9.49 




R 3 Bi 


7.2 


7.27 


7.66 


6.89 






7.11 


7.53 


7.61 


7.02 




V x Bi u 


7.60 


7.80 


8.07 


7.52 




B}B lg 


7.80 


7.87 


8.07 


7.59 




R 1 B 2s 


8.01 


8.15 


8.19 


7.66 




R!A 9 


8.29 


8.36 


8.52 


7.87 


C2H4 




8.62 


8.76 


8.80 


8.36 




V 3 Bi u 


4.36 


4.64 


4.99 


4.12 




R 3 B 3ll 


6.98 


7.43 


7.48 


6.92 




R 3 Bi 9 


7.79 


7.47 


7.82 


7.50 




R 3 B 29 


7.79 


8.06 


8.07 


7.56 




R 3 A 9 


8.15 


8.12 


8.11 


7.63 




V 1 A 2 


4.07 


3.63 


3.68 


3.88 




R 1 B 2 


7.11 


7.48 


7.92 


6.96 




R X B 2 


7.97 


8.13 


8.58 


7.66 




R : Ai 


8.14 


9.13 


9.47 


8.74 




R X A 2 


8.37 


8.30 


8.84 


7.84 


CH2O 


R!B 2 


8.88 


8.86 


9.21 


8.52 




V 3 A 2 


3.50 


3.02 


3.12 


3.21 




V 3 A X 


5.86 


5.70 


6.02 


5.29 
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FIG. 5. The lowest CT excitation of C 2 H 4 - • • C 2 F 4 dimer along the intermolecular distance R (in A). For all functionals, 
the excitation at 5 A is set to zero. 



Since cjM05-D is a parametrized functional, we test it ag ainst the trial simple model wM05s-D as well as two 
closely related functionals (M05-2X 19] and o;B97X-D [36]) on a separate independent test set of data, which 
includes further atomization energies, reaction energies, noncovalent interaction energies, equilibrium geometries, 
energy curve for homonuclear diatomic cation dissociations, frontier orbital energies and fundamental gaps. The 
three databases assessing frontier orbital energies and fundamental gaps are presented for the first time. For exci- 
tation energies, we calculate valence and Rydberg excitations, as well as a charge-transfer excited state. Compared 
to wM05s-D, noticeable difference in transferability for atomization energies largely decides our proposed model. 
cjM05-D consistently outperforms M05-2X (and performs comparably to cjB97X-D) on the test sets, and shows 
smaller SIE and better asymptotic behavior relative to both the M05-2X and wB97X-D. 
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FIG. 6. The lowest CT excitation of C 2 H 4 - • • C 2 F 4 dimer along the intermolecular distance R (in A). 
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SUPPLEMENTARY MATERIAL 



TABLE XIII: -e N (N) (in eV) of IP131 database. Most 
experimental reference values are collected from the NIST 
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.72 


14. 


.28 


BC1 3 (Borane, trichloro-) 


11, 


.64 


11. 


,18 


10. 


.58 


11. 


02 


AlClg (Aluminum trichloride) 


12 


.01 


11. 


.48 


10. 


.85 


11. 


.29 


CF 4 (Carbon tetrafluoride) 


16. 


.20 


15. 


23 


15 


.15 


14. 


.69 


CC1 4 (Carbon tetrachloride) 


11. 


69 


11. 


.15 


10. 


.58 


10. 


99 


OCS (Carbonyl sulfide) 


11. 


19 


10. 


.64 


9. 


.96 


10. 


53 


CS 2 (Carbon disulfide) 


10. 


.09 


9. 


.57 


8. 


.89 


9. 


.57 


CF 2 (Carbonic difluoride) 


13. 


.60 


12. 


.78 


12. 


.43 


12. 


.38 


SiF 4 (Silicon tetrafluoride) 


16. 


.40 


15. 


37 


15 


26 


14. 


.85 


N 2 (Nitrous oxide) 


12. 


.89 


12. 


00 


11. 


.51 


11. 


.75 


NF 3 (Nitrogen trifluoride) 


13. 


.60 


12. 


62 


12 


35 


12 


32 


PF 3 (Phosphorus trifluoride) 


12. 


.20 


10. 


.74 


10. 


.17 


10. 


69 


3 (Ozone) 


12. 


73 


12. 


.35 


11. 


.91 


11. 


.97 


F 2 (Difluorine monoxide) 


13. 


.26 


12. 


51 


12 


.40 


12 


.02 


CIF3 (Chlorine trifluoride) 


13 


.05 


12. 


19 


11. 


.86 


11. 


.86 


C 2 F 4 (Tetrafluoroethylene) 


10. 


.69 


9. 


.93 


9. 


,41 


9. 


.71 


CF 3 CN (Acetonitrile, trifluoro-) 


14. 


.30 


13. 


30 


12. 


,73 


13 


.00 


CH 3 CCH (Propyne) 


10. 


.37 


9 


.82 


9. 


.08 


9, 


.60 


CH 2 CCH 2 (Allene) 


10. 


.20 


9. 


.79 


9. 


.06 


9. 


60 


C 3 H 4 (Cyclopropene) 


9. 


.86 


9. 


,33 


8. 


.60 


9. 


.11 


C 3 H 6 (Cyclopropane) 


10. 


.54 


10. 


.47 


9. 


.74 


10. 


.28 


C 3 H g (Propane) 


11. 


51 


11, 


.21 


10. 


.66 


11. 


.10 


CH3CCCH3 (2-Butyne) 


9. 


.79 


9. 


19 


8. 


.46 


8. 


.98 


C 4 H 6 (Cyclobutene) 


9. 


.43 


9. 


.19 


8. 


.43 


8, 


.98 


CH 3 CH(CH 3 )CH3 (Isobutane) 


11. 


.13 


10. 


99 


10. 


.42 


10. 


.85 


C 6 H 6 (Benzene) 


9. 


.25 


9. 


.28 


8. 


.49 


9, 


06 


CH 2 F 2 (Methane, difluoro-) 


13. 


.27 


12. 


.27 


11. 


97 


11, 


.97 


CHF 3 (Methane, trifluoro-) 


15. 


50 


13. 


.52 


13. 


,27 


13. 
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CH 2 C1 2 (Methylene chloride) 

CHCI3 (Chloroform) 

CH 3 N0 2 (Methane, nitro-) 

CH 3 SiH 3 (Methyl silane) 

HCOOH (Formic acid) 

CH 3 CONH 2 (Acetamide) 

C 2 H 5 N (Aziridine) 

C 2 N 2 (Cyanogen) 

CH 3 NHCH 3 (Dimethylamine) 

CH 2 CO (Ketene) 

C 2 H 4 (Ethylene oxide) 

C 2 H 2 2 (Ethanedial) 

CH 3 CH 2 OH (Ethanol) 

CH 3 OCH 3 (Dimethyl ether) 

C 2 H 4 S (Thiirane) 

CH 3 SOCH 3 (Dimethyl sulfoxide) 

CH 2 CHF (Ethene, fluoro-) 

CH 3 CH 2 C1 (Ethyl chloride) 

CH 2 CHC1 (Ethene, chloro-) 

CH 3 COCl (Acetyl Chloride) 

CH 2 C1CH 2 CH 3 (Propane, 1-chloro-) 

N(CH 3 ) 3 (Trimethylamine) 

C 4 H 4 (Furan) 

C 4 H 5 N (Pyrrole) 

N0 2 (Nitrogen dioxide) 

SF 6 (Sulfur Hexafluoride) [87] 

CFC1 3 (Trichloromonofluoromethane) 

CF 3 C1 (Methane, chlorotrifluoro-) 

CF 3 Br (Bromotrifluoromethane) [87] 

HCCF (Fluoroacetylene) [87] 

HCCCN (Cyanoacetylene) [87] 



11.40 


10.85 


11.50 


10.85 


11.29 


11.12 


11.60 


11.23 


11.50 


10.69 


10.00 


9.66 


9.85 


9.36 


13.51 


12.78 


8.95 


8.60 


9.64 


9.19 


10.57 


10.25 


10.60 


10.12 


10.64 


10.09 


10.10 


9.66 


9.05 


8.57 


9.10 


8.76 


10.63 


9.90 


11.06 


10.44 


10.20 


9.60 


11.03 


10.77 


10.88 


10.42 


8.54 


8.27 


8.90 


8.70 


8.23 


8.13 


11.23 


10.72 


15.70 


14.96 


11.76 


11.23 


13.08 


12.27 


12.08 


11.29 


11.50 


10.53 


11.75 


11.07 



10.28 


10.72 


10.28 


10.72 


10.77 


10.72 


10.50 


11.12 


10.23 


10.36 


9.17 


9.36 


8.73 


9.11 


12.19 


12.54 


7.97 


8.38 


8.49 


9.00 


9.79 


9.90 


9.63 


9.87 


9.63 


9.79 


9.14 


9.36 


7.86 


8.46 


8.16 


8.57 


9.17 


9.66 


9.85 


10.28 


8.87 


9.41 


10.23 


10.53 


9.76 


10.25 


7.62 


8.05 


7.92 


8.51 


7.32 


7.92 


10.42 


10.36 


14.91 


14.42 


10.66 


11.07 


11.75 


12.05 


10.74 


11.12 


9.85 


10.28 


10.42 


10.85 
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C 4 N 2 (2-Butynedinitrile) [87] 


11 


,84 


11.56 


10.96 


11.34 


C 2 N 2 (Cyanogen) 


13 


.51 


12.78 


12.19 


12.54 


C 3 2 (Carbon suboxide) [87] 


10 


,80 


10.39 


9.82 


10.23 


FCN (Cyanogen fluoride) [87] 


13 


.65 


12.46 


11.89 


12.16 


C 4 H 2 (Diacetylene) 


10 


.30 


9.71 


9.00 


9.52 


H 2 CS (Thioformaldehyde) 


9 


.38 


8.73 


8.02 


8.62 


CHONH 2 (Formamide) [87] 


10 


.40 


9.90 


9.44 


9.57 


CH 2 CHCHO (Acrolein) [87] 


10 


.10 


9.90 


9.38 


9.57 


CH 2 CC1 2 (Ethene, 1,1-dichloro-) 


10 


.00 


9.57 


8.87 


9.41 


C 2 HF 3 (Trifluoroethylene) [87] 


10 


.62 


9.74 


9.14 


9.52 


CH 2 CF 2 (Ethene, 1,1-difluoro-) [87] 


10 


.70 


10.01 


9.33 


9.76 


CH 3 F (Methyl fluoride) 


13 


.04 


12.21 


11.86 


11.91 


CF 2 C1 2 (Difluorodichloromethane) 


12 


.24 


11.61 


11.07 


11.45 


SiF 2 (Silicon difluoride) 


11 


.08 


10.25 


9.52 


10.23 


MSE 






-0.79 


-1.34 


-1.02 


MAE 






0.79 


1.34 


1.02 


rms 






0.96 


1.43 


1.18 



TABLE XIV: -e7v+i(iV + 1) and -e N+ i(N) (in eV) of EA115 
database. Reference values are calculated by the CCSD(T) 
theory. 



molecule 


reference 


wM05-D 
-l(iV + l) -e N+ i(N) - ejV4 


M05-2X 

-i(N+l) -e N+1 (N) - ejV4 


wB97X-D 

-i(JV + l) -e N +l(N) 


H (Hydrogen atom) 


0.75 


0.84 


-0.52 


-0.08 


-0.53 


0.84 


-0.53 


He (Helium atom) 


-2.63 


-4.35 


-4.76 


-5.17 


-3.89 


-4.27 


-4.30 


Li (Lithium atom) 


0.62 


0.68 


-0.44 


0.00 


0.14 


0.71 


-0.19 


Be (Beryllium atom) 


-0.36 


-0.46 


-0.52 


-1.20 


0.35 


-0.41 


-0.30 


B (Boron atom) 


0.25 


0.19 


0.38 


-0.65 


1.28 


0.16 


0.57 


C (Carbon atom) 


1.25 


0.87 


1.69 


0.08 


2.42 


0.73 


1.99 


N (Nitrogen atom) 


-0.22 


-0.60 


0.49 


-1.39 


1.55 


-0.63 


0.52 


O (Oxygen atom) 


1.45 


0.76 


2.18 


0.16 


2.96 


0.54 


2.48 
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F (Fluorine atom) 


3.44 


2.28 


4.30 


1.90 


4.81 


1.93 


4.90 


Ne (Neon atom) 


-5.31 


-7.18 


-7.07 


-8.16 


-5.82 


-7.53 


-6.83 


Na (Sodium atom) 


0.54 


0.68 


-0.49 


0.05 


-0.03 


0.68 


-0.24 


Mg (Magnesium atom) 


-0.23 


-0.41 


-0.65 


-0.98 


0.11 


-0.30 


-0.38 


Al (Aluminum atom) 


0.45 


0.27 


0.22 


-0.54 


1.12 


0.33 


0.35 


Si (Silicon atom) 


1.42 


1.09 


1.41 


0.22 


2.31 


1.06 


1.55 


P (Phosphorus atom) 


0.74 


0.57 


0.98 


-0.24 


2.20 


0.63 


1.20 


S (Sulfur atom) 


2.10 


1.77 


2.53 


1.01 


3.64 


1.71 


2.77 


CI (Chlorine atom) 


3.69 


3.05 


4.27 


2.34 


5.22 


2.91 


4.57 


Ar (Argon atom) 


-2.81 


-3.18 


-3.78 


-4.03 


-2.53 


-3.48 


-3.40 


CH 3 (Methyl radical) 


-0.07 


-0.08 


-0.08 


-0.92 


1.09 


-0.16 


-0.05 


CH 4 (Methane) 


-0.62 


-1.01 


-1.52 


-1.58 


-0.92 


-0.92 


-1.17 


NH (Imidogen) 


0.33 


-0.03 


0.65 


-0.79 


1.66 


-0.14 


0.84 


NH 2 (Amino radical) 


0.74 


0.44 


0.90 


-0.33 


1.85 


0.27 


1.12 


NH 3 (Ammonia) 


-0.56 


-0.98 


-1.52 


-1.63 


-0.76 


-0.90 


-1.12 


OH (Hydroxyl radical) 


1.83 


1.22 


2.15 


0.63 


2.91 


0.95 


2.53 


H 2 (Water) 


-0.56 


-1.03 


-1.50 


-1.80 


-0.68 


-0.95 


-1.09 


HF (Hydrogen fluoride) 


-0.63 


-1.14 


-1.47 


-1.85 


-0.73 


-1.01 


-1.12 


SiH 3 (Silyl) 


0.93 


0.79 


1.14 


-0.03 


2.53 


0.82 


1.31 


SiH 4 (Silane) 


-1.11 


-1.09 


-1.69 


-1.60 


-0.95 


-1.06 


-1.28 


PH 3 (Phosphine) 


-1.21 


-0.90 


-1.50 


-1.44 


-0.71 


-0.90 


-1.09 


SH 2 (Hydrogen sulfide) 


-0.49 


-1.85 


-1.41 


-2.39 


-0.63 


-1.69 


-1.03 


HC1 (Hydrogen chloride) 


-0.52 


-0.92 


-1.36 


-1.63 


-0.54 


-0.84 


-0.98 


C 2 H 2 (Acetylene) 


-1.90 


-2.50 


-1.52 


-3.24 


-1.09 


-2.56 


-1.22 


C 2 H 4 (Ethylene) 


-1.86 


-2.07 


-1.74 


-2.86 


-0.92 


-2.12 


-1.33 


C 2 H 6 (Ethane) 


-0.62 


-1.03 


-1.55 


-1.55 


-0.90 


-0.95 


-1.17 


HCN (Hydrogen cyanide) 


-0.48 


-2.07 


-1.41 


-2.83 


-0.92 


-2.01 


-1.09 


CO (Carbon monoxide) 


-1.50 


-1.90 


-1.33 


-2.67 


-0.46 


-1.93 


-1.09 


HCO (Formyl radical) 


0.02 


-0.27 


0.44 


-1.01 


1.39 


-0.33 


0.57 


CH 2 (Formaldehyde) 


-0.55 


-1.31 


-0.63 


-2.09 


0.19 


-1.41 


-0.46 


CH 3 OH (Methyl alcohol) 


-0.55 


-0.95 


-1.44 


-1.50 


-0.79 


-0.84 


-1.09 
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N 2 (Nitrogen diatomic) 


-2.24 


-2.80 


-1.55 


-3.56 


-0.71 


-2.88 


-1.28 


N 2 H 4 (Hydrazine) 


-0.45 


-1.52 


-1.41 


-1.99 


-0.57 


-1.39 


-1.01 


NO (Nitric oxide) 


-0.42 


-0.92 


0.60 


-1.52 


1.44 


-1.12 


0.87 


2 (Oxygen diatomic) 


-0.08 


-0.87 


0.79 


-1.36 


1.60 


-1.09 


1.12 


H 2 2 (Hydrogen peroxide) 


-0.92 


-1.47 


-1.63 


-2.04 


-0.82 


-1.33 


-1.25 


F 2 (Fluorine diatomic) 


0.42 


-0.38 


1.69 


-0.57 


2.42 


-0.73 


2.15 


C0 2 (Carbon dioxide) 


-0.65 


-4.49 


-1.69 


-5.11 


-0.73 


-4.68 


-1.31 


P 2 (Phosphorus diatomic) 


0.48 


0.33 


0.73 


-0.46 


1.71 


0.35 


0.95 


S 2 (Sulfur diatomic) 


1.53 


1.12 


1.69 


0.44 


2.72 


1.12 


1.93 


Cl 2 (Chlorine diatomic) 


0.75 


0.38 


1.06 


-0.33 


1.99 


0.27 


1.33 


NaCl (Sodium Chloride) 


0.65 


0.52 


0.16 


-0.16 


0.79 


0.65 


0.46 


SiO (Silicon monoxide) 


0.03 


-0.33 


0.11 


-1.12 


1.09 


-0.24 


0.33 


CS (Carbon monosulfide) 


-0.09 


-0.30 


0.33 


-1.03 


1.25 


-0.33 


0.54 


CIO (Monochlorine monoxide) 


2.19 


1.60 


2.58 


1.06 


3.37 


1.41 


2.77 


C1F (Chlorine monofluoride) 


0.44 


-0.14 


0.92 


-0.79 


1.85 


-0.22 


1.28 


Si 2 H 6 (Disilane) 


-0.69 


-1.22 


-1.58 


-1.69 


-0.76 


-1.17 


-1.20 


CH3CI (Methyl chloride) 


-0.51 


-0.84 


-1.39 


-1.47 


-0.68 


-0.76 


-1.01 


CH 3 SH (Methanethiol) 


-0.50 


-0.87 


-1.41 


-1.44 


-0.63 


-0.79 


-1.03 


S0 2 (Sulfur dioxide) 


0.81 


0.41 


1.41 


-0.19 


2.42 


0.30 


1.66 


BF 3 (Borane, trifluoro-) 


-1.04 


-1.25 


-1.71 


-1.88 


-0.82 


-1.17 


-1.31 


BCI3 (Borane, trichloro-) 


-0.17 


-1.01 


-0.16 


-1.52 


0.76 


-0.98 


-0.03 


AICI3 (Aluminum trichloride) 


0.06 


-0.08 


-0.35 


-0.79 


0.73 


-0.11 


-0.19 


CF 4 (Carbon tetrafluoride) 


-1.33 


-1.96 


-2.48 


-2.53 


-1.52 


-1.82 


-2.04 


CC1 4 Carbon tetrachloride) 


-0.46 


-0.35 


-0.41 


-1.01 


0.57 


-0.35 


-0.14 


OCS (Carbonyl sulfide) 


-0.74 


-1.58 


-1.03 


-2.26 


-0.08 


-1.58 


-0.82 


CS 2 (Carbon disulfide) 


0.01 


-0.11 


0.16 


-0.79 


1.12 


-0.11 


0.35 


CF 2 (Carbonic difluoride) 


-2.37 


-2.50 


-1.55 


-3.40 


-0.65 


-2.42 


-1.36 


SiF 4 (Silicon tetrafluoride) 


-0.81 


-1.25 


-1.66 


-1.88 


-0.63 


-1.22 


-1.22 


N 2 (Nitrous oxide) 


-2.01 


-2.58 


-1.58 


-3.26 


-0.65 


-2.75 


-1.36 


NF 3 (Nitrogen trifluoride) 


-2.06 


-2.88 


-2.67 


-3.64 


-1.82 


-2.86 


-2.37 


PF 3 (Phosphorus trifluoride) 


-1.23 


-1.58 


-1.88 


-2.12 


-0.98 


-1.50 


-1.55 



3 (Ozone) 


1.93 


1.90 


F 2 (Difluorine monoxide) 


-0.31 


-1.01 


C1F 3 (Chlorine trifluoride) 


1.20 


0.63 


C 2 F 4 (Tetrafluoroethylene) 


-1.65 


-2.42 


CH 3 CCH (Propyne) 


-1.13 


-1.77 


CH 2 CCH 2 (Allene) 


-0.56 


-1.55 


C 3 H 4 (Cyclopropene) 


-1.82 


-2.01 


C 3 H 6 (Cyclopropane) 


-0.65 


-1.17 


CH 2 F 2 (Methane, difluoro-) 


-0.58 


-1.06 


CF 3 H (Methane, trifluoro-) 


-0.60 


-1.17 


CH 2 C1 2 (Methylene chloride) 


-0.49 


-0.82 


CHCI3 (Chloroform) 


-0.83 


-0.79 


CH 3 N0 2 (Methane, nitro-) 


-0.37 


-0.49 


CH 3 SiH 3 (Methyl silane) 


-0.53 


-0.92 


HCOOH (Formic acid) 


-0.57 


-2.15 


CH 3 CONH 2 (Acetamide) 


-0.31 


-1.60 


C 2 N 2 (Cyanogen) 


-0.19 


-0.33 


CH 2 CO (Ketene) 


-0.51 


-1.41 


C 2 H 4 (Ethylene oxide) 


-0.86 


-1.03 


C 2 H 2 2 (Ethanedial) 


0.69 


0.63 


CH 3 CH 2 OH (Ethanol) 


-0.53 


-0.92 


CH3OCH3 (Dimethyl ether) 


-0.58 


-1.01 


C 2 H 4 S (Thiirane) 


-0.78 


-1.25 


CH 2 CHF (Ethene, fluoro-) 


-0.88 


-2.18 


CH 3 CH 2 C1 (Ethyl chloride) 


-0.51 


-0.90 


CH 2 CHC1 (Ethene, chloro-) 


-1.11 


-1.60 


CH3CCIO (Acetyl Chloride) 


-0.85 


-1.12 


N0 2 (Nitrogen dioxide) 


1.44 


0.90 


CFCI3 (Trichloromonofluoromethane) 


-0.68 


-0.76 


CF3CI (Methane, chlorotrifluoro-) 


-1.06 


-1.63 


HCCF (Fluoroacetylene) 


-0.55 


-0.95 



3.24 


1.58 


4.22 


1.55 


3.48 


0.44 


-1.36 


1.22 


-1.17 


0.82 


1.60 


0.14 


2.56 


0.52 


1.90 


-2.07 


-2.94 


-1.03 


-2.34 


-1.66 


-1.41 


-2.23 


-0.82 


-1.60 


-1.06 


-1.60 


-2.04 


-0.92 


-1.41 


-1.22 


-1.66 


-2.80 


-0.92 


-2.07 


-1.31 


-1.71 


-1.63 


-1.03 


-1.09 


-1.33 


-1.52 


-1.63 


-0.95 


-0.90 


-1.20 


-1.58 


-1.85 


-1.01 


-1.01 


-1.25 


-1.31 


-1.50 


-0.49 


-0.73 


-0.98 


-1.01 


-1.52 


-0.05 


-0.73 


-0.76 


0.05 


-1.06 


0.98 


-0.65 


0.27 


-1.47 


-1.47 


-0.73 


-0.84 


-1.09 


-1.63 


-2.88 


-0.87 


-2.18 


-1.25 


-1.22 


-2.09 


-0.49 


-1.44 


-0.84 


0.52 


-1.03 


1.47 


-0.38 


0.68 


-1.03 


-1.96 


-0.16 


-1.28 


-0.87 


-1.60 


-1.55 


-0.95 


-0.95 


-1.22 


1.28 


-0.11 


2.18 


0.49 


1.41 


-1.44 


-1.47 


-0.76 


-0.82 


-1.09 


-1.52 


-1.47 


-0.84 


-0.90 


-1.14 


-1.60 


-1.69 


-0.87 


-1.14 


-1.20 


-1.66 


-2.91 


-1.01 


-2.23 


-1.28 


-1.47 


-1.44 


-0.73 


-0.79 


-1.09 


-1.41 


-2.39 


-0.52 


-1.66 


-1.20 


-0.79 


-1.80 


0.11 


-1.14 


-0.65 


1.99 


0.38 


3.02 


0.76 


2.18 


-0.65 


-1.39 


0.33 


-0.73 


-0.41 


-1.82 


-2.28 


-0.84 


-1.55 


-1.52 


-1.39 


-1.47 


-0.82 


-0.82 


-1.09 
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HCCCN (Cyanoacetylene) 


-0.36 


-0.76 


-0.19 


-1.52 


0.71 


-0.79 


-0.08 


C 4 N 2 (2-Butynedinitrile) 


0.68 


0.76 


1.25 


0.05 


2.20 


0.71 


1.36 


C 2 N 2 (Cyanogen) 


-0.19 


-0.33 


0.52 


-1.03 


1.47 


-0.38 


0.68 


C 3 2 (Carbon suboxide) 


-0.74 


-1.01 


-0.19 


-1.66 


0.76 


-1.12 


-0.05 


FCN (Cyanogen fluoride) 


-0.66 


-4.08 


-1.31 


-4.81 


-0.33 


-4.22 


-0.95 


C 4 N 2 (Diacetylene) 


-0.64 


-1.28 


-0.92 


-2.01 


-0.03 


-1.28 


-0.82 


H 2 CS (Thioformaldehyde) 


0.28 


0.11 


0.68 


-0.68 


1.55 


0.05 


0.82 


CHONH 2 (Formamide) 


-0.35 


-2.48 


-1.33 


-3.21 


-0.63 


-2.45 


-0.95 


CH 2 CC1 2 (Ethene, 1,1-dichloro-) 


-1.07 


-1.01 


-1.12 


-1.58 


-0.19 


-0.90 


-0.98 


C 2 HF 3 (Trifluoroethylene) 


-0.54 


-1.14 


-1.60 


-1.74 


-0.92 


-0.98 


-1.25 


CH 2 CF 2 (Ethene, 1,1-difluoro-) 


-1.03 


-2.26 


-1.66 


-2.96 


-0.95 


-2.31 


-1.28 


CH 3 F (Methyl fluoride) 


-0.58 


-0.98 


-1.50 


-1.55 


-0.90 


-0.87 


-1.14 


CF 2 C1 2 (Difluorodichloromethane) 


-0.90 


-1.31 


-1.20 


-1.90 


-0.22 


-1.28 


-0.92 


SiF 2 (Silicon difluoride) 


0.10 


-0.27 


0.22 


-1.06 


1.20 


-0.19 


0.44 


MSE 




-0.53 


-0.24 


-1.18 


0.60 


-0.54 


0.04 


MAE 




0.55 


0.62 


1.18 


0.79 


0.56 


0.55 


rms 




0.79 


0.73 


1.32 


0.97 


0.82 


0.65 



TABLE XV: Reference values of FG115 database calculated 
by the CCSD(T) theory, compared with HOMO-LUMO gap 
results by DFT methods. 



Molecule 


Reference 


wM05-D 


M05-2X 


wB97X-D 


H (Hydrogen atom) 


12.86 


11.62 


10.98 


11.57 


He (Helium atom) 


27.23 


25.57 


24.70 


24.56 


Li (Lithium atom) 


4.22 


5.49 


4.05 


5.41 


Be (Beryllium atom) 


9.66 


8.70 


6.94 


8.57 


B (Boron atom) 


7.99 


6.83 


5.11 


6.66 


C (Carbon atom) 


9.97 


8.05 


6.69 


7.59 


N (Nitrogen atom) 


14.74 


11.97 


10.47 


11.61 


O (Oxygen atom) 


12.14 


9.63 


8.54 


9.00 


F (Fluorine atom) 


13.98 


10.88 


10.28 


9.74 
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Ne (Neon atom) 


26.91 


25.84 


24.75 


24.89 


Na (Sodium atom) 


4.14 


5.41 


4.11 


5.11 


Mg (Magnesium atom) 


7.76 


7.81 


6.26 


7.45 


Al (Aluminum atom) 


5.53 


5.11 


3.32 


5.11 


Si (Silicon atom) 


6.73 


5.88 


4.11 


5.74 


P (Phosphorus atom) 


9.78 


8.35 


6.31 


8.05 


S (Sulfur atom) 


8.23 


6.80 


5.00 


6.53 


CI (Chlorine atom) 


9.30 


7.43 


5.90 


6.99 


Ar (Argon atom) 


18.65 


17.95 


16.21 


17.33 


CH 3 (Methyl radical) 


9.86 


8.87 


6.83 


8.54 


CH 4 (Methane) 


15.06 


14.69 


13.52 


14.14 


NH (Imidogen) ( 3 XT) 


13.17 


11.32 


9.85 


10.80 


NH 2 (Amino radical) 


11.34 


10.09 


8.62 


9.63 


NH 3 (Ammonia) 


11.54 


11.53 


10.17 


10.85 


OH (Hydroxyl radical) 


11.27 


9.49 


8.40 


8.73 


H 2 (Water) 


13.35 


13.00 


11.78 


12.19 


HF (Hydrogen fluoride) 


16.91 


15.88 


15.01 


15.01 


SiH 3 (Silyl) 


7.95 


7.04 


4.79 


6.85 


SiH 4 (Silane) 


14.03 


13.55 


12.13 


13.08 


PH 3 (Phosphine) 


11.82 


11.29 


9.74 


10.85 


H 2 S (Hydrogen sulfide) 


11.00 


10.99 


9.47 


10.50 


HC1 (Hydrogen chloride) 


13.36 


13.00 


11.56 


12.43 


C 2 H 2 (Acetylene) 


13.43 


12.13 


10.93 


11.59 


C 2 H 4 (Ethylene) 


12.57 


11.75 


10.09 


11.07 


C 2 H 6 (Ethane) 


13.41 


13.25 


12.05 


12.76 


HCN (Hydrogen cyanide) 


14.31 


14.04 


12.89 


13.44 


CO (Carbon monoxide) 


15.57 


14.31 


12.97 


13.90 


HCO (Formyl radical) 


9.56 


8.38 


6.91 


8.05 


CH 2 (Formaldehyde) 


11.56 


10.72 


9.41 


10.28 


CH 3 OH (Methyl alcohol) 


11.67 


11.67 


10.55 


11.02 


N 2 (Nitrogen diatomic) 


17.88 


16.13 


14.96 


15.56 
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N 2 H 4 (Hydrazine) 


10. 


29 


10 


.44 


9. 


.03 


9. 


.79 


NO (Nitric oxide) 


10. 


.11 


8. 


.11 


6. 


.91 


7. 


.48 


2 (Oxygen diatomic) ( 3 £ 9 ) 


12. 


52 


10. 


.69 


9. 


.66 


9. 


90 


H 2 2 (Hydrogen peroxide) 


12. 


.65 


12. 


.40 


11. 


23 


11. 


.61 


F 2 (Fluorine diatomic) 


15 


53 


12. 


.76 


11. 


.97 


11. 


.75 


C0 2 (Carbon dioxide) 


14. 


.58 


14. 


.66 


13. 


.27 


13 


.95 


P 2 (Phosphorus diatomic) 


10. 


.19 


9. 


.19 


7. 


.43 


8, 


.95 


S 2 (Sulfur diatomic) ( 3 S g ) 


7, 


.96 


7. 


.26 


5. 


.58 


6. 


.94 


Cl 2 (Chlorine diatomic) 


10 


93 


9. 


.71 


8. 


.19 


9. 


.30 


NaCl (Sodium Chloride) 


8. 


.64 


8. 


.43 


7. 


.15 


7. 


.94 


SiO (Silicon monoxide) 


11. 


.60 


10. 


.77 


9. 


.11 


10. 


.44 


CS (Carbon monosulfide) 


11. 


.58 


10. 


66 


9 


.17 


10. 


.36 


CIO (Monochlorine monoxide) 


8. 


.85 


7. 


53 


6. 


.34 


7. 


13 


C1F (Chlorine monofluoride) 


12 


.43 


10. 


.74 


9. 


.33 


10 


.17 


Si 2 H 6 (Disilane) 


11. 


33 


11. 


.72 


10 


.17 


11. 


32 


CH3CI (Methyl chloride) 


12. 


01 


12. 


.00 


10. 


66 


11. 


.42 


CH 3 SH (Methanethiol) 


10. 


.01 


10. 


.20 


8, 


.68 


9. 


.71 


S0 2 (Sulfur dioxide) 


11. 


.74 


10. 


36 


8. 


.87 


9. 


93 


BF 3 (Borane, trifluoro-) 


17. 


.22 


16. 


.54 


15 


53 


15 


.59 


BCI3 (Borane, trichloro-) 


12. 


.07 


11. 


.34 


9. 


.82 


11. 


.04 


AICI3 (Aluminum trichloride) 


12. 


.13 


11. 


83 


10. 


.12 


11. 


.48 


CF 4 (Carbon tetrafluoride) 


17 


.85 


17. 


.71 


16 


.67 


16. 


73 


CC1 4 (Carbon tetrachloride) 


11. 


.97 


11. 


.56 


10. 


01 


11. 


.12 


OCS (Carbonyl sulfide) 


12. 


13 


11. 


.67 


10. 


.04 


11. 


.34 


CS 2 (Carbon disulfide) 


10. 


19 


9. 


.41 


7. 


.78 


9. 


.22 


CF 2 (Carbonic difluoride) 


16. 


.08 


14. 


33 


13. 


.08 


13. 


.74 


SiF 4 (Silicon tetrafluoride) 


16. 


.95 


17. 


03 


15 


.88 


16. 


.08 


N 2 (Nitrous oxide) 


15. 


.01 


13. 


.57 


12. 


.16 


13. 


.11 


NF 3 (Nitrogen trifluoride) 


15. 


.76 


15. 


.29 


14. 


.17 


14. 


69 


PF 3 (Phosphorus trifluoride) 


13. 


00 


12 


.62 


11. 


.15 


12. 


.24 


3 (Ozone) 


11. 


.06 


9. 


.11 


7. 


.70 


8. 


49 
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F 2 (Difluorine monoxide) 


13 


.82 


12 


.08 


11. 


.18 


11. 


.21 


ClFg (Chlorine trifluoride) 


11. 


.79 


10 


.58 


9. 


.30 


9 


.96 


C 2 F 4 (Tetrafluoroethylene) 


12 


.45 


12 


,00 


10. 


.44 


11 


.37 


CH 3 CCH (Propyne) 


11. 


.69 


11 


.23 


9. 


.90 


10 


,66 


CH 2 CCH 2 (Allene) 


10 


.83 


11. 


.40 


9. 


.98 


10 


.83 


C 3 H 4 (Cyclopropene) 


11 


,87 


10 


.99 


9. 


.52 


10 


,42 


C 3 H 6 (Cyclopropane) 


11 


,64 


12 


.19 


10 


,77 


11. 


,61 


CH 2 F 2 (Methane, difluoro-) 


14 


.15 


13 


.79 


12. 


92 


13 


,16 


CF 3 H (Methane, trifluoro-) 


15. 


.44 


15 


,10 


14. 


.28 


14 


,47 


CH 2 C1 2 (Methylene chloride) 


12 


.18 


12 


,16 


10. 


.77 


11. 


.70 


CHCI3 (Chloroform) 


12 


.38 


11 


.86 


10. 


.34 


11. 


,48 


CH 3 N0 2 (Methane, nitro-) 


11. 


.94 


11. 


.07 


9. 


.79 


10 


,44 


CH 3 SiH 3 (Methyl silane) 


12 


.35 


12 


.70 


11. 


23 


12 


,21 


HCOOH (Formic acid) 


11 


.98 


12 


.32 


11, 


.10 


11 


,61 


CH 3 CONH 2 (Acetamide) 


10 


.05 


10 


.88 


9. 


.66 


10 


.20 


C 2 N 2 (Cyanogen) 


13 


.90 


12 


.27 


10. 


.72 


11. 


.86 


CH 2 CO (Ketene) 


10 


.32 


10 


.23 


8. 


.65 


9. 


,87 


C 2 H 4 (Ethylene oxide) 


11 


,68 


11 


.86 


10. 


.74 


11. 


.12 


C 2 H 2 2 (Ethanedial) 


10 


.04 


8 


.84 


7, 


.45 


8. 


,46 


CH 3 CH 2 OH (Ethanol) 


11. 


.38 


11 


.53 


10. 


39 


10 


,88 


CH3OCH3 (Dimethyl ether) 


10 


,79 


11. 


.18 


9 


.98 


10 


,50 


C 2 H 4 S (Thiirane) 


9 


.93 


10 


,17 


8, 


.73 


9 


.66 


CH 2 CHF (Ethene, fluoro-) 


11. 


.55 


11. 


,56 


10. 


.17 


10 


.93 


CH 3 CH 2 C1 (Ethyl chloride) 


11. 


.74 


11 


.91 


10. 


58 


11. 


.37 


CH 2 CHC1 (Ethene, chloro-) 


11. 


.35 


11 


,02 


9. 


.38 


10 


,61 


CH3COCI (Acetyl Chloride) 


11 


,97 


11. 


.56 


10. 


.12 


11. 


,18 


N0 2 (Nitrogen dioxide) 


9 


.79 


8 


.73 


7. 


.40 


8. 


,19 


CFCI3 (Trichloromonofluoromethane) 


12 


,61 


11 


.89 


10. 


.34 


11. 


,48 


CF3CI (Methane, chlorotrifluoro-) 


14 


.27 


14 


.09 


12. 


59 


13 


.57 


HCCF (Fluoroacetylene) 


12 


.04 


11 


.91 


10. 


66 


11. 


.37 


HCCCN (Cyanoacetylene) 


12 


.20 


11 


,26 


9. 


.71 


10 


.93 
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C 4 N 2 (2-Butynedinitrile) 


11 


.52 


10.31 


8.76 


9.98 


C 2 N 2 (Cyanogen) 


13 


.90 


12.27 


10.72 


11.86 


C 3 2 (Carbon suboxide) 


11 


.64 


10.58 


9.06 


10.28 


FCN (Cyanogen fluoride) 


14 


.33 


13.76 


12.21 


13.11 


C 4 H 2 (Diacetylene) 


11 


,00 


10.64 


9.03 


10.34 


H 2 CS (Thioformaldehyde) 


9 


.18 


8.05 


6.47 


7.81 


CHONH 2 (Formamide) 


10 


.81 


11.23 


10.06 


10.53 


CH 2 CC1 2 (Ethene, 1,1-dichloro-) 


11 


.17 


10.69 


9.06 


10.39 


C 2 HF 3 (Trifluoroethylene) 


11 


.11 


11.34 


10.06 


10.77 


CH 2 CF 2 (Ethene, 1,1-difluoro-) 


11 


.81 


11.67 


10.28 


11.04 


CH 3 F (Methyl fluoride) 


14 


,09 


13.71 


12.76 


13.06 


CF 2 C1 2 (Difluorodichloromethane) 


13 


.33 


12.81 


11.29 


12.38 


SiF 2 (Silicon difluoride) 


11 


.04 


10.04 


8.32 


9.79 


MSE 






-0.70 


-2.08 


-1.21 


MAE 






0.84 


2.08 


1.25 


rms 






1.07 


2.24 


1.48 
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TABLE XVI: Experimental molecular geometries (in 
angstrom) for IP131, EA115 and FG115 databases collected 



from the NIST database 



861 ] . The integers following the 



molecule are in order the spin multiplicities of neutral species, 
cation and anion. The molecules with only the spin multi- 
plicity of the neutral species shown are in neither EA115 nor 
FG115 database. The spin multiplicities of neutral species 
are experimental values specified in the NIST database 86]. 



Each spin multiplicity of cations and anions is decided by 
CCSD(T)/aug-cc-pVTZ calculations, and corresponds to the 
lowest energy in between various spin multiplicities. 



H (Hydrogen atom) 2 


Null 1 


H 






He (Helium atom) 1 


2 


2 


He 






Li (Lithium atom) 2 


1 


1 


Li 






Be (Beryllium atom) 


1 


2 2 


Be 






B (Boron atom) 2 1 


3 




B 






C (Carbon atom) 3 


2 


4 


C 






N (Nitrogen atom) 4 


3 


3 


N 






(Oxygen atom) 3 


4 


2 









F (Fluorine atom) 2 


3 


1 


F 







Ne (Neon atom) 12 2 
Ne 
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Na (Sodium atom) 2 11 
Na 

Mg (Magnesium atom) 12 2 

Mg 

Al (Aluminum atom) 2 13 
Al 

Si (Silicon atom) 3 2 4 
Si 

P (Phosphorus atom) 4 3 3 
P 

S (Sulfur atom) 3 4 2 
S 

CI (Chlorine atom) 2 3 1 
CI 



Ar (Argon atom) 12 2 
Ar 



CH 3 (Methyl radical) 2 1 


1 






CI 


0.0000 


0.0000 


0.0000 


H2 


1.0790 


0.0000 


0.0000 


H3 


-0.5395 


-0.9344 


0.0000 


H4 


-0.5395 


0.9344 


0.0000 


CH 4 (Methane) 1 2 2 








CI 


0.0000 


0.0000 


0.0000 


H2 


0.6276 


0.6276 


0.6276 


H3 


0.6276 


-0.6276 


-0.6276 


H4 


-0.6276 


0.6276 


-0.6276 


H5 


-0.6276 


-0.6276 


0.6276 


NH (Imidogen) ( 3 XT) 3 2 


2 






Nl 


0.0000 


0.0000 


0.0000 


H2 


0.0000 


0.0000 


1.0362 
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NH 2 (Amino radical) 2 3 1 



Nl 0.0000 0.0000 0.0000 

H2 0.0000 0.8036 0.6347 

H3 0.0000 -0.8036 0.6347 
NH 3 (Ammonia) 12 2 

Nl 0.0000 0.0000 0.0000 

H2 0.0000 -0.9377 -0.3816 

H3 0.8121 0.4689 -0.3816 

H4 -0.8121 0.4689 -0.3816 
OH (Hydroxyl radical) 2 3 1 

01 0.0000 0.0000 0.0000 

H2 0.0000 0.0000 0.9697 

H 2 (Water) 1 2 2 

01 0.0000 0.0000 0.1173 

H2 0.0000 0.7572 -0.4692 

H3 0.0000 -0.7572 -0.4692 
HF (Hydrogen fluoride) 12 2 

Fl 0.0000 0.0000 0.0000 

H2 0.0000 0.0000 0.9168 
SiH 3 (Silyl) 2 1 1 

Sil 0.0000 0.0000 0.0819 

H2 0.0000 1.3928 -0.3820 

H3 1.2062 -0.6964 -0.3820 

H4 -1.2062 -0.6964 -0.3820 

SiH 4 (Silane) 1 2 2 

Sil 0.0000 0.0000 0.0000 

H2 0.8544 0.8544 0.8544 

H3 -0.8544 -0.8544 0.8544 

H4 -0.8544 0.8544 -0.8544 
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H5 0.8544 -0.8544 -0.8544 

PH 3 (Phosphine) 1 2 2 

PI 0.0000 0.0000 0.0000 

H2 0.0000 -1.1932 -0.7717 

H3 1.0333 0.5966 -0.7717 

H4 -1.0333 0.5966 -0.7717 
H 2 S (Hydrogen sulfide) 1 2 2 

SI 0.0000 0.0000 0.0000 

H2 0.0000 0.9569 0.9208 

H3 0.0000 -0.9569 0.9208 
HC1 (Hydrogen chloride) 12 2 

Cll 0.0000 0.0000 0.0000 

H2 0.0000 0.0000 1.2746 

C 2 H 2 (Acetylene) 12 2 

CI 0.0000 0.0000 0.6013 

C2 0.0000 0.0000 -0.6013 

H3 0.0000 0.0000 1.6644 

H4 0.0000 0.0000 -1.6644 

C 2 H 4 (Ethylene) 1 2 2 

CI 0.0000 0.0000 0.6695 

C2 0.0000 0.0000 -0.6695 

H3 0.0000 0.9289 1.2321 

H4 0.0000 -0.9289 1.2321 

H5 0.0000 0.9289 -1.2321 

H6 0.0000 -0.9289 -1.2321 
C 2 H 6 (Ethane) 12 2 

CI 0.0000 0.0000 0.7680 

C2 0.0000 0.0000 -0.7680 

H3 -1.0192 0.0000 1.1573 
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H4 0.5096 


0.8826 


1.1573 


H5 0.5096 


-0.8826 


1.1573 


H6 1.0192 


0.0000 


-1.1573 


H7 -0.5096 


-0.8826 


-1.1573 


H8 -0.5096 


0.8826 


-1.1573 


HCN (Hydrogen cyanide) 12 2 






CI 0.0000 


0.0000 


0.0000 


H2 0.0000 


0.0000 


1.0640 


N3 0.0000 


0.0000 


-1.1560 


CO (Carbon monoxide) 12 2 






CI 0.0000 


0.0000 


0.0000 


02 0.0000 


0.0000 


1.1283 


HCO (Formyl radical) 2 1 1 






CI 0.0000 


0.0000 


0.0000 


H2 1.0800 


0.0000 


0.0000 


03 -0.5899 


1.0427 


0.0000 


H 2 CO (Formaldehyde) 1 2 2 






01 0.0000 


0.0000 


1.2050 


C2 0.0000 


0.0000 


0.0000 


H3 0.0000 


0.9429 


-0.5876 


H4 0.0000 


-0.9429 


-0.5876 


CH 3 OH (Methyl alcohol) 1 2 2 






CI -0.0503 


0.6685 


0.0000 


02 -0.0503 


-0.7585 


0.0000 


H3 -1.0807 


1.0417 


0.0000 


H4 0.4650 


1.0417 


0.8924 


H5 0.4650 


1.0417 


-0.8924 


H6 0.8544 


-1.0677 


0.0000 



N 2 (Nitrogen diatomic) 12 2 
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Nl 0.0000 


0.0000 


0.5488 


N2 0.0000 


0.0000 


-0.5488 


N 2 H 4 (Hydrazine) 12 2 






Nl 0.0000 


0.7230 


-0.1123 


N2 0.0000 


-0.7230 


-0.1123 


H3 -0.4470 


1.0031 


0.7562 


H4 0.4470 


-1.0031 


0.7562 


H5 0.9663 


1.0031 


0.0301 


H6 -0.9663 


-1.0031 


0.0301 


NO (Nitric oxide) 2 1 3 






01 0.0000 


0.0000 


0.0000 


N2 0.0000 


0.0000 


1.1508 


2 (Oxygen diatomic) ( 3 S ff ) 3 2 2 






01 0.0000 


0.0000 


0.0000 


02 0.0000 


0.0000 


1.2075 


H 2 2 (Hydrogen peroxide) 12 2 






1 0.0000 


0.7375 


-0.0528 


02 0.0000 


-0.7375 


-0.0528 


H3 0.8190 


0.8170 


0.4220 


H4 -0.8190 


-0.8170 


0.4220 


F 2 (Fluorine diatomic) 12 2 






Fl 0.0000 


0.0000 


0.0000 


F2 0.0000 


0.0000 


1.4119 


C0 2 (Carbon dioxide) 12 2 






CI 0.0000 


0.0000 


0.0000 


02 0.0000 


0.0000 


1.1621 


03 0.0000 


0.0000 


-1.1621 



P 2 (Phosphorus diatomic) 12 2 
PI 0.0000 



0.0000 



0.0000 
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P2 0.0000 0.0000 1.8934 

S 2 (Sulfur diatomic) ( 3 S 9 ) 3 2 2 

51 0.0000 0.0000 0.0000 

52 0.0000 0.0000 1.8892 
Cl 2 (Chlorine diatomic) 12 2 

Cll 0.0000 0.0000 0.0000 

C12 0.0000 0.0000 1.9879 
NaCl (Sodium Chloride) 1 2 2 

Nal 0.0000 0.0000 0.0000 

C12 0.0000 0.0000 2.3608 

510 (Silicon monoxide) 12 2 

511 0.0000 0.0000 0.0000 
02 0.0000 0.0000 1.5097 
CS (Carbon monosulfide) 12 2 

CI 0.0000 0.0000 0.0000 

S2 0.0000 0.0000 1.5349 

CIO (Monochlorine monoxide) 2 3 1 

01 0.0000 0.0000 0.0000 

C12 0.0000 0.0000 1.5696 

C1F (Chlorine monofluoride) 12 2 

Fl 0.0000 0.0000 0.0000 

C12 0.0000 0.0000 1.6283 
Si 2 H 6 (Disilane) 12 2 

511 0.0000 0.0000 1.1600 

512 0.0000 0.0000 -1.1600 
H3 0.0000 1.3865 1.6483 
H4 -1.2008 -0.6933 1.6483 
H5 1.2008 -0.6933 1.6483 
H6 0.0000 -1.3865 -1.6483 
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H7 -1.2008 


0.6933 


-1.6483 


H8 1.2008 


0.6933 


-1.6483 


CH 3 C1 (Methyl chloride) 1 2 2 
CI 0.0000 


0.0000 


0.0000 


C12 0.0000 


0.0000 


1.7810 


H3 1.0424 


0.0000 


-0.3901 


H4 -0.5212 


0.9027 


-0.3901 


H5 -0.5212 


-0.9027 


-0.3901 


CH 3 SH (Methanethiol) 1 2 2 

CI -0.8500 


-0.0344 


-0.2000 


S2 0.9000 


-0.5125 


-0.1219 


H3 1.4219 


0.5781 


0.4250 


H4 -0.9406 


0.8688 


-0.8219 


H5 -1.4219 


-0.8688 


-0.6469 


H6 -1.2031 


0.1656 


0.8219 


S0 2 (Sulfur dioxide) 1 2 2 

SI 0.0000 


0.0000 


0.0000 


02 0.0000 


1.2371 


0.7215 


03 0.0000 


-1.2371 


0.7215 


BF 3 (Borane, trifluoro-) 12 2 
Bl 0.0000 


0.0000 


0.0000 


F2 0.0000 


1.3070 


0.0000 


F3 1.1319 


-0.6535 


0.0000 


F4 -1.1319 


-0.6535 


0.0000 


BC1 3 (Borane, trichloro-) 12 2 






Bl 0.0000 


0.0000 


0.0000 


C12 0.0000 


1.7421 


0.0000 


C13 1.5087 


-0.8711 


0.0000 


C14 -1.5087 


-0.8711 


0.0000 
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A1C1 3 (Aluminum trichloride) 12 2 






All 0.0000 


0.0000 


0.0000 


C12 0.0000 


2.0600 


0.0000 


C13 1.7840 


-1.0300 


0.0000 


C14 -1.7840 


-1.0300 


0.0000 


CF 4 (Carbon tetrafluoride) 12 2 






CI 0.0000 


0.0000 


0.0000 


F2 0.7593 


0.7593 


0.7593 


F3 -0.7593 


-0.7593 


0.7593 


F4 -0.7593 


0.7593 


-0.7593 


F5 0.7593 


-0.7593 


-0.7593 


CC1 4 (Carbon tetrachloride) 12 2 






CI 0.0000 


0.0000 


0.0000 


C12 1.0202 


1.0202 


1.0202 


C13 -1.0202 


-1.0202 


1.0202 


C14 -1.0202 


1.0202 


-1.0202 


C15 1.0202 


-1.0202 


-1.0202 


OCS (Carbonyl sulfide) 1 2 2 






CI 0.0000 


0.0000 


0.0000 


02 0.0000 


0.0000 


1.1600 


S3 0.0000 


0.0000 


-1.5600 


CS 2 (Carbon disulfide) 1 2 2 






CI 0.0000 


0.0000 


0.0000 


S2 0.0000 


0.0000 


1.5540 


S3 0.0000 


0.0000 


-1.5540 


CF 2 (Carbonic difluoride) 1 2 2 






01 0.0000 


0.0000 


1.3143 


C2 0.0000 


0.0000 


0.1403 


F3 0.0000 


1.0614 


-0.6309 
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F4 0.0000 


-1.0614 


-0.6309 


SiF 4 (Silicon tetrafluoride) 12 2 
Sil 0.0000 


0.0000 


0.0000 


F2 0.8972 


0.8972 


0.8972 


F3 -0.8972 


-0.8972 


0.8972 


F4 -0.8972 


0.8972 


-0.8972 


F5 0.8972 


-0.8972 


-0.8972 


N 2 (Nitrous oxide) 1 2 2 

Nl 0.0000 


0.0000 


-1.1998 


N2 0.0000 


0.0000 


-0.0716 


03 0.0000 


0.0000 


1.1126 


NF 3 (Nitrogen trifluoride) 12 2 
Nl 0.0000 


0.0000 


0.4731 


F2 0.0000 


1.2279 


-0.1226 


F3 1.0634 


-0.6140 


-0.1226 


F4 -1.0634 


-0.6140 


-0.1226 


PF 3 (Phosphorus trifluoride) 12 2 
PI 0.0000 


0.0000 


0.4602 


F2 0.0000 


1.3578 


-0.2557 


F3 1.1759 


-0.6789 


-0.2557 


F4 -1.1759 


-0.6789 


-0.2557 


3 (Ozone) 1 2 2 

01 0.0000 


0.0000 


0.0000 


02 0.0000 


1.0885 


0.6697 


03 0.0000 


-1.0885 


0.6697 


F 2 (Difluorine monoxide) 12 2 
1 0.0000 


0.0000 


0.6074 


F2 0.0000 


1.1063 


-0.2700 


F3 0.0000 


-1.1063 


-0.2700 



54 



C1F 3 (Chlorine trifluoride) 


1 2 2 






Cll 


0.0000 


0.0000 


0.3572 


F2 


0.0000 


0.0000 


-1.2408 


F3 


0.0000 


1.6964 


0.2831 


F4 


0.0000 


-1.6964 


0.2831 


C 2 F 4 (Tetrafluoroethylene) 


1 2 2 






CI 


0.0000 


0.0000 


0.6555 


C2 


0.0000 


0.0000 


-0.6555 


F3 


0.0000 


1.0961 


1.3893 


F4 


0.0000 


-1.0961 


1.3893 


F5 


0.0000 


-1.0961 


-1.3893 


F6 


0.0000 


1.0961 


-1.3893 


CF 3 CN (Acetonitrile, trifluoro-) 1 






CI 


0.0000 


0.0000 


-0.3406 


C2 


0.0000 


0.0000 


1.1518 


N3 


0.0000 


0.0000 


2.3054 


F4 


0.0000 


1.2540 


-0.7780 


F5 


1.0860 


-0.6270 


-0.7780 


F6 


-1.0860 


-0.6270 


-0.7780 


CH 3 CCH (Propyne) 1 2 


2 






CI 


0.0000 


0.0000 


-1.2455 


C2 


0.0000 


0.0000 


0.2135 


C3 


0.0000 


0.0000 


1.4195 


H4 


0.0000 


0.0000 


2.4755 


H5 


0.0000 


1.0465 


-1.6003 


H6 


0.9063 


-0.5232 


-1.6003 


H7 


-0.9063 


-0.5232 


-1.6003 


CH 2 CCH 2 (Allene) 1 2 


2 






CI 


0.0000 


0.0000 


0.0000 
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C2 


0.0000 


0.0000 


1.3080 


C3 


0.0000 


0.0000 


-1.3080 


H4 


0.0000 


0.9327 


1.8662 


H5 


0.0000 


-0.9327 


1.8662 


H6 


0.9327 


0.0000 


-1.8662 


H7 


-0.9327 


0.0000 


-1.8662 


C 3 H 4 (Cyclopropene) 


1 2 2 






CI 


0.0000 


0.0000 


0.8628 


C2 


0.0000 


0.6476 


-0.5001 


C3 


0.0000 


-0.6476 


-0.5001 


H4 


0.0000 


1.5745 


-1.0386 


H5 


0.0000 


-1.5745 


-1.0386 


H6 


0.9154 


0.0000 


1.4509 


H7 


-0.9154 


0.0000 


1.4509 


C 3 H 6 (Cyclopropane) 


12 2 






CI 


0.0000 


0.8666 


0.0000 


C2 


0.7505 


-0.4333 


0.0000 


C3 


-0.7505 


-0.4333 


0.0000 


H4 


0.0000 


1.4525 


0.9108 


H5 


1.2579 


-0.7262 


0.9108 


H6 


-1.2579 


-0.7262 


0.9108 


H7 


0.0000 


1.4525 


-0.9108 


H8 


1.2579 


-0.7262 


-0.9108 


H9 


-1.2579 


-0.7262 


-0.9108 


C 3 H g (Propane) 1 








CI 


0.0000 


0.5863 


-0.0000 


C2 


-1.2681 


-0.2626 


0.0000 


C3 


1.2681 


-0.2626 


-0.0000 


H4 


0.0000 


1.2449 


0.8760 



56 



H5 


-0.0003 


1.2453 


-0.8758 


H6 


-2.1576 


0.3742 


0.0000 


H7 


2.1576 


0.3743 


-0.0000 


H8 


-1.3271 


-0.9014 


0.8800 


H9 


-1.3271 


-0.9014 


-0.8800 


H10 


1.3271 


-0.9014 


-0.8800 


Hll 


1.3272 


-0.9014 


0.8800 


CH3CCCH3 (2-Butyne) 1 








CI 


0.0000 


0.0000 


0.6070 


C2 


0.0000 


0.0000 


-0.6070 


C3 


0.0000 


0.0000 


2.0750 


C4 


0.0000 


0.0000 


-2.0750 


H5 


0.0000 


1.0440 


2.4695 


H6 


-0.9041 


-0.5220 


2.4695 


H7 


0.9041 


-0.5220 


2.4695 


H8 


0.0000 


1.0440 


-2.4695 


H9 


0.9041 


-0.5220 


-2.4695 


HIO 


-0.9041 


-0.5220 


-2.4695 


C 4 H 6 (Cyclobutene) 1 








CI 


0.0000 


0.6710 


0.8107 


C2 


0.0000 


-0.6710 


0.8107 


C3 


0.0000 


0.7821 


-0.7023 


C4 


0.0000 


-0.7821 


-0.7023 


H5 


0.0000 


1.4165 


1.5962 


H6 


0.0000 


-1.4165 


1.5962 


H7 


0.8986 


1.2425 


-1.1233 


H8 


-0.8986 


-1.2425 


-1.1233 


H9 


-0.8986 


1.2425 


-1.1233 


HIO 


0.8986 


-1.2425 


-1.1233 
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CH 3 CH(CH 3 )CH 3 (Isobutane) 


1 






CI 


0.0000 


0.0000 


0.3650 


H2 


0.0000 


0.0000 


1.4730 


C3 


0.0000 


1.4528 


-0.0987 


C4 


1.2582 


-0.7264 


-0.0987 


C5 


-1.2582 


-0.7264 


-0.0987 


H6 


0.0000 


1.4867 


-1.1931 


H7 


1.2875 


-0.7433 


-1.1931 


H8 


-1.2875 


-0.7433 


-1.1931 


H9 


0.8941 


1.9575 


0.2821 


H10 


-0.8941 


1.9575 


0.2821 


Hll 


1.2482 


-1.7530 


0.2821 


H12 


2.1422 


-0.2045 


0.2821 


H13 


-2.1422 


-0.2045 


0.2821 


H14 


-1.2482 


-1.7530 


0.2821 


C 6 H 6 (Benzene) 1 








CI 


0.0000 


1.3970 


0.0000 


C2 


1.2098 


0.6985 


0.0000 


C3 


1.2098 


-0.6985 


0.0000 


C4 


0.0000 


-1.3970 


0.0000 


C5 


-1.2098 


-0.6985 


0.0000 


C6 


-1.2098 


0.6985 


0.0000 


H7 


0.0000 


2.4810 


0.0000 


H8 


2.1486 


1.2405 


0.0000 


H9 


2.1486 


-1.2405 


0.0000 


H10 


0.0000 


-2.4810 


0.0000 


Hll 


-2.1486 


-1.2405 


0.0000 


H12 


-2.1486 


1.2405 


0.0000 



CH 2 F 2 (Methane, dinuoro-) 12 2 
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CI 0.0000 


0.0000 


0.5003 


H2 -0.9029 


0.0000 


1.1002 


H3 0.9029 


0.0000 


1.1002 


F4 0.0000 


1.0962 


-0.2890 


F5 0.0000 


-1.0962 


-0.2890 


CHF 3 (Methane, trifluoro-) 12 2 
CI 0.0000 


0.0000 


0.3086 


H2 0.0000 


0.0000 


1.4066 


F3 0.0000 


1.2609 


-0.1207 


F4 1.0920 


-0.6305 


-0.1207 


F5 -1.0920 


-0.6305 


-0.1207 


CH 2 C1 2 (Methylene chloride) 1 2 2 
CI 0.0000 


0.0000 


0.7761 


H2 -0.8854 


0.0000 


1.3734 


H3 0.8854 


0.0000 


1.3734 


C14 0.0000 


1.4675 


-0.2178 


C15 0.0000 


-1.4675 


-0.2178 


CHC1 3 (Chloroform) 1 2 2 

ci 0.0000 


0.0000 


0.5231 


H2 0.0000 


0.0000 


1.5961 


C13 0.0000 


1.6562 


-0.0928 


C14 I.4343 


-0.8281 


-0.0928 


C15 -1.4343 


-0.8281 


-0.0928 


CH 3 N0 2 (Methane, nitro-) 12 2 
CI 1.4008 


0.0000 


0.0000 


N2 -0.0878 


0.0000 


0.0000 


H3 1.7215 


-1.0392 


0.0000 


H4 1.7215 


0.5196 


0.9000 


H5 1.7215 


0.5196 


-0.9000 



59 



06 


-0.6498 


1.0874 


0.0000 


07 


-0.6498 


-1.0874 


0.0000 


CH 3 SiH 3 (Methyl silane) 1 


2 2 






CI 


0.0000 


0.0000 


-1.2367 


Si2 


0.0000 


0.0000 


0.6319 


H3 


0.0000 


-1.0237 


-1.6272 


H4 


-0.8866 


0.5119 


-1.6272 


H5 


0.8866 


0.5119 


-1.6272 


H6 


0.0000 


1.3893 


1.1514 


H7 


-1.2031 


-0.6946 


1.1514 


H8 


1.2031 


-0.6946 


1.1514 


HC00H (Formic acid) 1 2 


2 






CI 


0.0000 


0.4199 


0.0000 


02 


-1.0543 


-0.4121 


0.0000 


03 


1.1506 


0.0721 


0.0000 


H4 


-0.0799 


1.5140 


0.0000 


H5 


-0.6905 


-1.3134 


0.0000 


CH 3 CONH 2 (Acetamide) 1 


2 2 






CI 


-1.3674 


-0.3302 


0.0013 


C2 


0.0720 


0.1552 


-0.0019 


N3 


1.0259 


-0.8416 


-0.0326 


04 


0.3726 


1.3376 


0.0057 


H5 


-2.0656 


0.5484 


-0.0612 


H6 


-1.5426 


-1.0080 


-0.8780 


H7 


-1.5719 


-0.9088 


0.9430 


H8 


2.0050 


-0.5675 


0.0710 


H9 


0.7853 


-1.8244 


0.1113 


C 2 H 5 N (Aziridine) 1 








Nl 


-0.0037 


0.8559 


-0.1934 
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H2 




-0.0140 


1.4040 


0.6620 


C3 




-0.7380 


-0.4043 


0.0261 


C4 




0.7425 


-0.3952 


0.0380 


H5 




-1.2757 


-0.4486 


0.9652 


H6 




-1.2332 


-0.8641 


-0.8216 


H7 




1.2820 


-0.7584 


-0.8280 


H8 




1.2404 


-0.5270 


0.9918 


C 2 N 2 (Cyanogen) 1 2 
CI 


2 


0.0000 


0.0000 


0.6950 


C2 




0.0000 


0.0000 


-0.6950 


N3 




0.0000 


0.0000 


1.8450 


N4 




0.0000 


0.0000 


-1.8450 


CH3NHCH3 (Dimethylamine) 
Nl 


1 

0.0934 


-0.5520 


0.0000 


H2 




-0.6454 


-1.2544 


0.0000 


C3 




-0.0220 


0.2554 


1.2139 


C4 




-0.0220 


0.2554 


-1.2139 


H5 




-0.9221 


0.8833 


1.2437 


H6 




-0.9221 


0.8833 


-1.2437 


H7 




0.8440 


0.9287 


1.2609 


H8 




0.8440 


0.9287 


-1.2609 


H9 




-0.0138 


-0.3909 


2.0843 


H10 




-0.0138 


-0.3909 


-2.0843 


CH 2 CO (Ketene) 1 2 
CI 


2 


0.0000 


0.0000 


0.0000 


C2 




0.0000 


0.0000 


1.3150 


03 




0.0000 


0.0000 


2.4750 


H4 




0.0000 


0.9451 


-0.5206 


H5 




0.0000 


-0.9451 


-0.5206 
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C 2 H 4 (Ethylene oxide) 1 


2 2 






01 


0.0000 


0.0000 


0.8517 


C2 


0.0000 


0.7297 


-0.3725 


C3 


0.0000 


-0.7297 


-0.3725 


H4 


-0.9230 


1.2565 


-0.5859 


H5 


0.9230 


1.2565 


-0.5859 


H6 


0.9230 


-1.2565 


-0.5859 


H7 


-0.9230 


-1.2565 


-0.5859 


C 2 H 2 2 (Ethanedial) 1 2 


2 






CI 


0.0000 


0.7630 


0.0000 


C2 


0.0000 


-0.7630 


0.0000 


H3 


1.0481 


1.1907 


0.0000 


H4 


-1.0481 


-1.1907 


0.0000 


05 


-1.0367 


1.3908 


0.0000 


06 


1.0367 


-1.3908 


0.0000 


CH 3 CH 2 OH (Ethanol) 1 2 


2 






CI 


1.1879 


-0.3829 


0.0000 


C2 


0.0000 


0.5526 


0.0000 


03 


-1.1867 


-0.2472 


0.0000 


H4 


-1.9237 


0.3850 


0.0000 


H5 


2.0985 


0.2306 


0.0000 


H6 


1.1184 


-1.0093 


0.8869 


H7 


1.1184 


-1.0093 


-0.8869 


H8 


-0.0227 


1.1812 


0.8852 


H9 


-0.0227 


1.1812 


-0.8852 


CH 3 OCH 3 (Dimethyl ether) 


1 2 2 






01 


0.0000 


0.0000 


0.5952 


C2 


0.0000 


1.1669 


-0.1963 


C3 


0.0000 


-1.1669 


-0.1963 
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H4 


0.0000 


2.0489 


0.4542 


H5 


0.0000 


-2.0489 


0.4542 


H6 


0.8950 


1.1787 


-0.8287 


H7 


-0.8950 


1.1787 


-0.8287 


H8 


-0.8950 


-1.1787 


-0.8287 


H9 


0.8950 


-1.1787 


-0.8287 


C 2 H 4 S (Thiirane) 1 2 


2 






SI 


0.0000 


0.0000 


0.8622 


C2 


0.0000 


0.7421 


-0.7942 


C3 


0.0000 


-0.7421 


-0.7942 


H4 


-0.9174 


1.2493 


-1.0661 


H5 


0.9174 


1.2493 


-1.0661 


H6 


0.9174 


-1.2493 


-1.0661 


H7 


-0.9174 


-1.2493 


-1.0661 


CH3SOCH3 (Dimethyl sulfoxide) 1 






SI 


0.0000 


0.1432 


0.4202 


02 


0.0000 


1.4024 


-0.3667 


C3 


1.3425 


-0.8664 


-0.2227 


C4 


-1.3425 


-0.8664 


-0.2227 


H5 


2.2553 


-0.4311 


0.0751 


H6 


-2.2553 


-0.4311 


0.0751 


H7 


1.3057 


-1.8579 


0.2349 


H8 


1.2255 


-0.9383 


-1.3113 


H9 


-1.3057 


-1.8579 


0.2349 


H10 


-1.2255 


-0.9383 


-1.3113 


CH 2 CHF (Ethene, fluoro- 


) 1 2 2 






Cl 


0.0000 


0.4476 


0.0000 


C2 


1.1877 


-0.1487 


0.0000 


F3 


-1.1356 


-0.2769 


0.0000 
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H4 


-0.2349 


1.5038 


0.0000 


H5 


1.2321 


-1.2348 


0.0000 


H6 


2.0966 


0.4290 


0.0000 


CH 3 CH 2 C1 (Ethyl chloride) 


1 2 2 






CI 


1.5949 


-0.3563 


-0.0000 


C2 


0.4757 


0.6568 


-0.0004 


H3 


2.5527 


0.1648 


0.0000 


H4 


1.5351 


-0.9919 


-0.8828 


H5 


1.5347 


-0.9917 


0.8828 


C16 


-1.1206 


-0.1505 


-0.0005 


H7 


0.5089 


1.2945 


-0.8790 


H8 


0.4949 


1.2796 


0.8893 


CH 2 CHC1 (Ethene, chloro-) 


1 2 2 






CI 


0.0000 


0.0000 


1.3320 


C2 


0.0000 


0.0000 


0.0000 


C13 


0.0000 


1.4589 


2.2543 


H4 


0.0000 


-0.9058 


1.9384 


H5 


0.0000 


0.9249 


-0.5557 


H6 


0.0000 


-0.9391 


-0.5313 


CHgCOCl (Acetyl Chloride) 


1 2 2 






CI 


0.0000 


0.5272 


0.0000 


C2 


1.4961 


0.6994 


0.0000 


03 


-0.8349 


1.3710 


0.0000 


C14 


-0.4665 


-1.2092 


0.0000 


H5 


1.7591 


1.7726 


0.0000 


H6 


1.9367 


0.2285 


0.8973 


H7 


1.9367 


0.2285 


-0.8973 


CH 2 C1CH 2 CH 3 (Propane, 1-chloro-) 1 






Cll 


1.7375 


0.1388 


0.0000 



C2 


0.0000 


0.5937 


0.0000 


H3 


-0.1639 


1.2199 


0.9054 


H4 


-0.1639 


1.2199 


-0.9054 


C5 


-0.8957 


-0.6405 


0.0000 


H6 


-0.6638 


-1.2607 


0.8947 


H7 


-0.6638 


-1.2607 


-0.8947 


C8 


-2.3710 


-0.2543 


0.0000 


H9 


-2.6338 


0.3458 


0.8998 


H10 


-2.6338 


0.3458 


-0.8998 


Hll 


-3.0133 


-1.1633 


0.0000 


N(CH 3 ) 3 (Trimethylamine) 


1 






Nl 


0.0000 


-0.0000 


-0.3852 


C2 


-0.9758 


-0.9757 


0.0634 


C3 


1.3330 


-0.3571 


0.0634 


C4 


-0.3572 


1.3328 


0.0634 


H5 


-1.0253 


-1.0252 


1.1492 


H6 


1.4006 


-0.3752 


1.1492 


H7 


-0.3753 


1.4004 


1.1492 


H8 


-1.9833 


-0.7091 


-0.3156 


H9 


-0.7093 


-1.9833 


-0.3156 


H10 


1.6060 


-1.3629 


-0.3156 


Hll 


2.0722 


0.3776 


-0.3156 


H12 


0.3774 


2.0722 


-0.3156 


H13 


-1.3631 


1.6057 


-0.3156 


C 4 H 4 (Furan) 1 








01 


0.0000 


0.0000 


1.1626 


C2 


0.0000 


1.0920 


0.3487 


C3 


0.0000 


-1.0920 


0.3487 


C4 


0.0000 


0.7169 


-0.9596 
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C5 


0.0000 


-0.7169 


-0.9596 


H6 


0.0000 


2.0473 


0.8439 


H7 


0.0000 


-2.0473 


0.8439 


H8 


0.0000 


1.3509 


-1.8290 


H9 


0.0000 


-1.3509 


-1.8290 


C 4 H 5 N (Pyrrole) 1 








Nl 


0.0000 


0.0000 


1.1218 


H2 


0.0000 


0.0000 


2.1178 


C3 


0.0000 


1.1209 


0.3341 


C4 


0.0000 


-1.1209 


0.3341 


C5 


0.0000 


0.7076 


-0.9847 


C6 


0.0000 


-0.7076 


-0.9847 


H7 


0.0000 


2.1084 


0.7614 


H8 


0.0000 


-2.1084 


0.7614 


H9 


0.0000 


1.3566 


-1.8429 


H10 


0.0000 


-1.3566 


-1.8429 


N0 2 (Nitrogen dioxide) 2 


1 1 






Nl 


0.0000 


0.0000 


0.0000 


02 


0.0000 


1.0989 


0.4653 


03 


0.0000 


-1.0989 


0.4653 


SF 6 (Sulfur Hexafluoride) 


1 






SI 


0.0000 


0.0000 


0.0000 


F2 


0.0000 


0.0000 


1.5607 


F3 


0.0000 


1.5607 


0.0000 


F4 


1.5607 


0.0000 


0.0000 


F5 


0.0000 


-1.5607 


0.0000 


F6 


-1.5607 


0.0000 


0.0000 


F7 


0.0000 


0.0000 


-1.5607 



CFClg (Trichloromonofluoromethane) 12 2 
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CI 0.0000 


0.0000 


0.2472 


F2 0.0000 


0.0000 


1.5922 


C13 0.0000 


1.6732 


-0.3101 


C14 1.4491 


-0.8366 


-0.3101 


C15 -1.4491 


-0.8366 


-0.3101 


CF 3 C1 (Methane, chlorotrifluoro-) 12 2 






ci 0.0000 


0.0000 


-0.3471 


C12 0.0000 


0.0000 


1.4049 


F3 0.0000 


1.2425 


-0.8074 


F4 1.0760 


-0.6212 


-0.8074 


F5 -1.0760 


-0.6212 


-0.8074 


CF 3 Br (Bromotrifluoromethane) 1 






ci 0.0000 


0.0000 


-0.8088 


Br2 0.0000 


0.0000 


1.1146 


F3 0.0000 


1.2455 


-1.2651 


F4 1.0787 


-0.6228 


-1.2651 


F5 -1.0787 


-0.6228 


-1.2651 


HCCF (Fluoroacetylene) 12 2 






ci 0.0000 


0.0000 


-0.0942 


C2 0.0000 


0.0000 


-1.2922 


F3 0.0000 


0.0000 


1.1848 


H4 0.0000 


0.0000 


-2.3452 


HCCCN (Cyanoacetylene) 1 2 2 






Nl 0.0000 


0.0000 


1.9018 


C2 0.0000 


0.0000 


0.7413 


C3 0.0000 


0.0000 


-0.6351 


C4 0.0000 


0.0000 


-1.8409 


H5 0.0000 


0.0000 


-2.9033 



C 4 N 2 (2-Butynedinitrile) 12 2 
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CI 


0.0000 


0.0000 


2.5100 


C2 


0.0000 


0.0000 


3.7000 


C3 


0.0000 


0.0000 


1.1400 


C4 


0.0000 


0.0000 


5.0700 


N5 


0.0000 


0.0000 


0.0000 


N6 


0.0000 


0.0000 


6.2100 


C 2 N 2 (Cyanogen) 12 2 








CI 


0.0000 


0.0000 


0.6950 


C2 


0.0000 


0.0000 


-0.6950 


N3 


0.0000 


0.0000 


1.8450 


N4 


0.0000 


0.0000 


-1.8450 


C 3 2 (Carbon suboxide) 1 


2 2 






CI 


0.0000 


0.0000 


0.0308 


C2 


0.0000 


1.2509 


0.0123 


C3 


0.0000 


-1.2509 


0.0123 


04 


0.0000 


2.3964 


-0.0207 


05 


0.0000 


-2.3964 


-0.0207 


FCN (Cyanogen fluoride) 1 


2 2 






CI 


0.0000 


0.0000 


0.0000 


F2 


0.0000 


0.0000 


-1.2620 


N3 


0.0000 


0.0000 


1.1590 


C 4 H 2 (Diacetylene) 12 2 








CI 


0.0000 


0.0000 


0.6890 


C2 


0.0000 


0.0000 


-0.6890 


C3 


0.0000 


0.0000 


1.8940 


C4 


0.0000 


0.0000 


-1.8940 


H5 


0.0000 


0.0000 


2.9520 


H6 


0.0000 


0.0000 


-2.9520 



H 2 CS (Thioformaldehyde) 1 2 2 
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SI 0.0000 0.0000 0.5846 

C2 0.0000 0.0000 -1.0262 

H3 0.0000 0.9244 -1.5980 

H4 0.0000 -0.9244 -1.5980 
CHONH 2 (Formamide) 1 2 2 

CI 0.0000 0.4165 0.0000 

02 1.1942 0.2217 0.0000 

N3 -0.9373 -0.5551 0.0000 

H4 -0.4299 1.4182 0.0000 

H5 -0.6608 -1.5171 0.0000 

H6 -1.9020 -0.2881 0.0000 
CH 2 CHCHO (Acrolein) 1 

CI -0.1496 -0.7423 0.0000 

C2 0.0000 0.7200 0.0000 

C3 1.2291 1.2662 0.0000 

04 -1.2313 -1.3044 0.0000 

H5 0.8007 -1.3120 0.0000 

H6 -0.9066 1.3142 0.0000 

H7 1.3716 2.3428 0.0000 

H8 2.1072 0.6273 0.0000 
CH 2 CC1 2 (Ethene, 1,1-dichloro) 12 2 

CI 0.0000 0.0000 1.7363 

C2 0.0000 0.0000 0.4123 

H3 0.0000 0.9266 2.2713 

H4 0.0000 -0.9266 2.2713 

C15 0.0000 1.4382 -0.5128 

C16 0.0000 -1.4382 -0.5128 
C 2 HF 3 (Trifluoroethylene) 12 2 

CI 0.0000 0.4436 0.0000 
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C2 -0.7139 -0.6916 0.0000 

F3 1.3153 0.4857 0.0000 

F4 -0.5506 1.6388 0.0000 

F5 -0.0873 -1.8784 0.0000 

H6 -1.8133 -0.7269 0.0000 
CH 2 CF 2 (Ethene, 1,1-difluoro) 1 2 2 

CI 0.0000 0.0000 1.3851 

C2 0.0000 0.0000 0.0701 

H3 0.0000 0.9419 1.9094 

H4 0.0000 -0.9419 1.9094 

F5 0.0000 1.0777 -0.6972 

F6 0.0000 -1.0777 -0.6972 
CH 3 F (Methyl fluoride) 1 2 2 

CI 0.0000 0.0000 -0.6289 

F2 0.0000 0.0000 0.7541 

H3 0.0000 1.0201 -1.0043 

H4 0.8835 -0.5101 -1.0043 

H5 -0.8835 -0.5101 -1.0043 
CF 2 C1 2 (Dinuorodichloromethane) 12 2 

CI 0.0000 0.0000 0.3171 

F2 0.0000 1.0758 1.1243 

F3 0.0000 -1.0758 1.1243 

C14 1.4505 0.0000 -0.6512 

C15 -1.4505 0.0000 -0.6512 
SiF 2 (Silicon difluoride) 1 2 2 

Sil 0.0000 0.0000 0.5703 

F2 0.0000 1.2249 -0.4436 

F3 0.0000 -1.2249 -0.4436 
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TABLE XVII. Atomic energies (in kcal/mol) from the H atom to the Ar atom 



3- 



Atom 


Reference 


wM05-D 


M05-2X 


wB97X-D 


H 


-313.76 


-311.52 


-313.31 


-315.51 


He 


-1822.12 


-1822.11 


-1823.25 


-1824.57 


Li 


-4692.56 


-4693.72 


-4694.37 


-4699.80 


Be 


-9203.92 


-9198.97 


-9204.50 


-9202.60 


B 


-15470.58 


-15466.17 


-15471.98 


-15468.00 


C 


-23748.12 


-23746.17 


-23752.06 


-23745.12 


N 


-34255.28 


-34254.32 


-34260.32 


-34251.37 





-47105.50 


-47107.44 


-47113.72 


-47103.64 


F 


-62584.04 


-62587.59 


-62593.80 


-62583.06 


Ne 


-80909.66 


-80910.31 


-80917.57 


-80906.81 


Na 


-101816.42 


-101814.09 


-101820.46 


-101818.06 


Mg 


-125535.30 


-125538.45 


-125549.12 


-125538.97 


Al 


-152074.59 


-152077.17 


-152084.44 


-152078.77 


Si 


-181575.73 


-181577.53 


-181581.20 


-181577.86 


P 


-214143.51 


-214141.04 


-214140.76 


-214140.27 


S 


-249818.09 


-249817.74 


-249815.58 


-249817.21 


(HI 


988747 ^7 


988747 D4 


98874^ 1 ^ 


98874R ^5 


Ar 


-331036.74 


-331035.35 


-331028.07 


-331036.25 


MSE 




0.37 


-3.01 


-0.05 


MAE 




2.02 


5.10 


2.51 


rms 




2.42 


6.31 


2.96 
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TABLE XVIII. Non-hydrogen transfer barrier heights (in kcal/mol) of the NHTBH38/04 set [63]. 



Reactions 




AE ref 


ojM05-D 


M05-2X 


ojB97X-D 


Heavy-atom transfer reactions 












H + N 2 0->OH + N 2 


yf 


18.14 


16.52 


18.82 


17.47 




\' r 


83.22 


80.54 


85.45 


77.75 


H + FH— >HF + H 


V f 


42.18 


38.34 


41.16 


40.56 




V 


42.18 


38.34 


41.16 


40.56 


H + C1H->HC1 + H 


yS 


18.00 


17.83 


19.64 


18.19 




V 


18.00 


17.83 


19.64 


18.19 


H + FCH 3 ->HF + CH 3 


yt 


30.38 


28.76 


31.70 


30.02 




v 


57.02 


54.26 


55.28 


54.47 


H + F 2 ->HF + F 


t rf 

v> 


2.27 


-0.07 


2.90 


-0.60 




v r 


106.18 


106.50 


114.64 


103.96 


CH 3 + FCl-s-CHgF + CI 


V' 


7.43 


4.22 


4.04 


2.76 




yr 


60.17 


57.63 


59.94 


56.55 


Nucleophilic substitution reactions 












F~ + CH 3 F->FCH 3 + F" 


V 1 


-0.34 


-1.79 


-1.93 


-1.37 




yr 


-0.34 


-1.79 


-1.93 


-1.37 


F~ • • ■ CH 3 F^-FCH 3 - ■ ■ F~ 


V 1 


13.38 


12.98 


14.01 


12.85 




y 


13.38 


12.98 


14.01 


12.85 


Cl" + CH 3 Cl->-ClCH 3 + Cl" 


yf 


3.10 


2.29 


0.75 


3.70 




V 


3.10 


2.29 


0.75 


3.70 


Cl~ ■ ■ ■ CH 3 Cl-s>ClCH 3 - ■ ■ Cl~ 


yf 


13.61 


12.85 


12.01 


14.30 




y 


13.61 


12.85 


12.01 


14.30 


F~ + CH 3 C1->FCH 3 + CI" 


V f 


-12.54 


-14.38 


-15.89 


-13.58 




V 


20.11 


21.65 


21.86 


21.26 


F~ • • • CH 3 C1->FCH 3 ■ ■ ■ Cl~ 




2.89 


2.41 


2.13 


3.07 




yr 


29.62 


31.11 


31.97 


30.50 


OH~ + CH 3 F-*-HOCH 3 + F" 


yi 


-2.78 


-3.96 


-4.22 


-3.31 




v 


17.33 


16.05 


15.54 


17.90 


OH~ ■ ■ ■ CH 3 F-*-HOCH 3 ■ ■ ■ F~ 


t rf 


10.96 


10.32 


11.14 


10.57 




\r 


47.20 


49.33 


48.84 


49.20 


Unimolecular and association reactions 












H + N 2 ->HN 2 




14.69 


11.94 


14.65 


12.32 




yr 


10.72 


12.01 


13.07 


13.40 


H + CO->HCO 


yi 


3.17 


2.91 


4.96 


3.37 




yr 


22.68 


24.38 


25.23 


26.20 


H + C2H4 — ^-CH^CH^ 


yi 


1.72 


2.46 


3.68 


3.03 




yr 


41.75 


44.44 


44.88 


45.37 


CHg + C2H4 — ^CH^CH^CH^ 




6.85 


4.13 


4.09 


4.57 




yr 


32.97 


33.17 


32.57 


33.85 


HCN -> HNC 


vf 


48.16 


44.78 


46.45 


46.42 




v~ 


33.11 


31.56 


33.10 


33.27 


MSE 






-0.94 


0.13 


-0.45 


MAE 






1.57 


1.75 


1.51 


rms 






1.89 


2.24 


2.00 
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TABLE XIX. Hydros 


$en transfer barrier heights (in 


kcal/mol) of the HTBH38/04 set Q 


r~i 

,y. 

7 ■— 'J 


Reactions 




AE ref 


ojM05-D 


M05-2X 


ojB97X-D 


H + HC1->H 2 + CI 


yf 


5.7 


3.61 


4.82 


4.08 




yr 


8.7 


4.68 


5.94 


4.58 


OH + H 2 ->H + H 2 


vf 


5.1 


1.94 


4.61 


2.22 




yr 


21.2 


18.69 


21.90 


18.77 


CH 3 + H 2 ->H + CH 4 


yf 


12.1 


9.14 


11.72 


9.06 




yr 


15.3 


13.35 


16.07 


13.54 


OH + CH 4 -s-CH 3 +H 2 


vf 


6.7 


2.53 


4.69 


3.04 




yr 


19.6 


15.06 


17.64 


15.11 


H + H 2 ->H 2 + H 


vf 


9.6 


8.83 


12.66 


9.21 




yr 


9.6 


8.83 


12.66 


9.21 


OH + NHg->H 2 + NH 2 


vf 


3.2 


-0.14 


2.67 


-0.02 




yr 


12.7 


9.73 


12.41 


9.69 


HC1 + CH 3 ->C1 + CH 4 


vf 


1.7 


-1.60 


-1.30 


-1.86 




yr 


7.9 


3.69 


4.17 


3.11 


OH + C 2 H 6 ^H 2 + C 2 H 5 


vf 


3.4 


-0.28 


2.19 


0.17 




yr 


19.9 


16.88 


19.24 


16.57 


F + H 2 ->HF + H 


vf 


1.8 


-2.22 


0.10 


-3.70 




yr 


33.4 


30.78 


33.40 


29.01 


O + CH 4 ^OH + CH 3 


vf 


13.7 


8.85 


10.62 


9.20 




yr 


8.1 


4.39 


5.50 


4.67 


H + PH 3 -i>PH 2 + H 2 


vf 


3.1 


2.21 


4.25 


3.28 




yr 


23.2 


23.45 


25.62 


23.46 


H + HO^H 2 + O 


vf 


10.7 


7.90 


10.38 


8.84 




yr 


13.1 


8.15 


11.16 


8.89 


H + H 2 S— >H 2 + HS 


vf 


3.5 


3.24 


5.42 


4.05 




yr 


17.3 


16.10 


18.21 


16.31 


O + HCl^OH + CI 


yf 


9.8 


4.95 


7.01 


5.47 




yr 


10.4 


5.78 


7.35 


5.91 


NH 2 + CH 3 ^CH 4 + NH 


vf 


8.0 


5.41 


7.10 


5.65 




yr 


22.4 


18.17 


20.33 


18.41 


NH 2 +C 2 H 5 ^C 2 H 6 +NH 


vf 


7.5 


7.09 


8.56 


7.26 




yr 


18.3 


15.22 


17.68 


15.69 


C 2 H 6 + NH 2 ^NH 3 + C 2 H 5 


vf 


10.4 


7.97 


10.11 


8.60 




yr 


17.4 


15.26 


17.42 


15.29 


NH 2 + CH 4 ^CH 3 + NH 3 


vf 


14.5 


10.81 


12.71 


11.28 




yr 


17.8 


13.47 


15.91 


13.65 


s— trans cis—C 5 H g —^s— trans cis— C 5 H g 


vf 


38.4 


37.23 


39.06 


39.07 




yr 


38.4 


37.23 


39.06 


39.07 


MSE 






-2.82 


-0.65 


-2.57 


MAE 






2.83 


1.51 


2.70 


rms 






3.16 


1.83 


3.10 
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TABLE XX. Interaction energies (in kcal/mol) for the S22 set with new reference values. [70| The counterpoise corrections 



are used to reduce the basis set superposition errors. Monomer deformation energies are not included. 



Complex [Symmetry] 


AE re f 


wM05-D 


M05-2X 


ojB97X-D 


Hydrogen bonded complexes 










(NH 3 ) 2 [C 2h ] 


-3.133 


-3.00 


-3.12 


-3.08 


(H 2 0) 2 [C s ] 


-4.989 


-4.89 


-5.08 


-4.97 


Formic acid dimer [C2/1] 


-18.753 


-18.98 


-19.41 


-19.28 


Formamide dimer [0%^] 


-16.062 


-15.65 


-15.81 


-16.17 


Uracil dimer [C2h] 


-20.641 


-19.90 


-19.64 


-20.43 


2-pyridoxine- 2-aminopyridine [C 1 ] 


-16.934 


-16.40 


-15.36 


-17.06 


Adenine-thymine WC [Ci] 


-16.660 


-15.87 


-14.84 


-16.46 


MSE 




0.35 


0.56 


-0.04 


MAE 




0.42 


0.77 


0.18 


Dispersion complexes 










(CH 4 ) 2 [D 3d ] 


-0.527 


-0.65 


-0.51 


-0.57 


(C 2 H 4 ) 2 [D M ] 


-1.472 


-1.70 


-1.42 


-1.79 


Benzene-CH 4 [C 3 ] 


-1.448 


-1.68 


-1.13 


-1.70 


Benzene dimer [C211] 


-2.654 


-3.55 


-1.38 


-3.16 


Pyrazine dimer [C s ] 


-4.255 


-4.40 


-2.95 


-4.27 


Uracil dimer [C2] 


-9.805 


-9.83 


-8.32 


-9.79 


Indolebenzene [Ci] 


-4.524 


-5.47 


-2.57 


-5.10 


Adenine-thymine stack [Ci] 


-11.730 


-12.00 


-9.40 


-11.86 


MSE 




-0.36 


1.09 


-0.23 


MAE 




0.36 


1.09 


0.23 


Mixed complexes 










Ethene-ethine [C2« ] 


-1.496 


-1.52 


-1.44 


-1.63 


Benzene-H20 [C s ] 


-3.275 


-3.60 


-3.57 


-3.50 


Benzene-NH3 [C s ] 


-2.312 


-2.54 


-2.23 


-2.53 


Benzene-HCN [C s ] 


-4.541 


-4.76 


-4.90 


-4.79 


Benzene dimer [C2 V ] 


-2.717 


-2.88 


-1.97 


-2.89 


Indole benzene T-shape [Ci] 


-5.627 


-5.85 


-4.76 


-5.64 


Phenol dimer [Ci] 


-7.097 


-6.89 


-6.03 


-6.99 


MSE 




-0.14 


0.31 


-0.13 


MAE 




0.20 


0.50 


0.16 


MSE 




-0.06 


0.67 


-0.14 


MAE 




0.33 


0.80 


0.19 


rms 




0.42 


1.06 


0.25 
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TABLE XXI. Comparison of errors of different functionals for the reaction energies (in kcal/mol) of the 30 chemical reactions 
in the NHTBH38/04 and HTBH38/04 database QQl- 



Reactions 


AE ref 


u;M05-D 


M05-2X 


u;B97X-D 


H + N 2 O^OH + N 2 


-65.08 


-64.03 


-66.63 


-60.28 


H + FCH 3 ^HF + CH 3 


-26.64 


-25.50 


-23.58 


-24.44 


H + F 2 ->HF + F 


-103.91 


-106.58 


-111.74 


-104.56 


CH 3 + FC1^CH 3 F + C1 


-52.74 


-53.41 


-55.90 


-53.79 


F" + CH 3 C1->FCH 3 + Cl- 


-32.65 


-36.03 


-37.75 


-34.84 


F" ■ ■ ■ CH 3 C1->FCH 3 • • • Cl" 


-26.73 


-28.71 


-29.84 


-27.43 


OH- + CH 3 F->HOCH 3 + F" 


-20.11 


-20.02 


-19.76 


-21.21 


OH- ■ ■ ■ CH 3 F-*-HOCH 3 - ■ ■ F" 


-36.24 


-39.01 


-37.70 


-38.63 


H + N 2 ->HN 2 


3.97 


-0.07 


1.58 


-1.09 


H + CO -> HCO 


-19.51 


-21.47 


-20.27 


-22.83 


H + C2H4 — ^CH^CH^ 


-40.03 


-41.98 


-41.20 


-42.35 


CH^ + C2H4 — ^-CHgCH^CH^ 


-26.12 


-29.04 


-28.48 


-29.28 


HON -> HNC 


15.05 


13.22 


13.35 


13.16 


H + HC1 -> H 2 + CI 


-3.0 


-1.08 


-1.12 


-0.49 


OH + H 2 -s>H + H 2 


-16.1 


-16.75 


-17.29 


-16.55 


CH 3 + H 2 ->H + CH 4 


-3.2 


-4.22 


-4.35 


-4.48 


OH + CH 4 -j-CH 3 + H 2 


-12.9 


-12.53 


-12.94 


-12.07 


OH + NH 3 ->H 2 + NH 2 


-9.5 


-9.87 


-9.75 


-9.71 


HC1 + CH 3 -s>Cl + CH 4 


-6.2 


-5.29 


-5.47 


-4.98 


OH + C 2 H 6 ^H 2 + C 2 H 5 


-16.5 


-17.16 


-17.06 


-16.40 


F + H 2 -> HF + H 


-31.6 


-33.00 


-33.29 


-32.71 


O + CH 4 -^OH + CH 3 


5.6 


4.46 


5.13 


4.53 


H + PH 3 ->PH 2 + H 2 


-20.1 


-21.25 


-21.37 


-20.18 


H + HO -v H 2 + O 


-2.4 


-0.25 


-0.78 


-0.05 


H + H 2 S->H 2 + HS 


-13.8 


-12.87 


-12.79 


-12.25 


O + HC1 -> OH + CI 


-0.6 


-0.83 


-0.34 


-0.45 


NH 2 + CH 3 ->CH 4 + NH 


-14.4 


-12.76 


-13.23 


-12.76 


NH 2 +C 2 H 5 ->C 2 H 6 + NH 


-10.8 


-8.13 


-9.11 


-8.43 


C 2 H 6 + NH 2 ^NH 3 +C 2 H 5 


-7.0 


-7.29 


-7.31 


-6.69 


NH 2 + CH 4 ^CH 3 + NH 3 


-3.3 


-2.66 


-3.20 


-2.37 


MSE 




-0.58 


-0.86 


-0.24 


MAE 




1.49 


1.65 


1.63 


rms 




1.79 


2.29 


2.06 
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TABLE XXI 
from Ref. 



92] 



Binding energies (in kcal/mol) of several sets of noncovalent interactions. The first three sets are taken 
with monomer deformation energies taken into considerations. The last three sets are taken from Ref. 64(] 



without considering monomer deformation energies. The counter-point corrections are applied for all the cases. 





Complex 


AB rB f 


qjM05-D 


M05-2X 


wB97X-D 


Charge-transfer complexes 










C 2 H 4 - ■ ■ F 2 


1.06 


0.79 


1.00 


0.78 


NH 3 ..-F 2 


1.81 


1.41 


1.44 


1.50 


C 2 H 2 - ■ ■ C1F 


3.81 


3.82 


4.51 


3.66 


HCN- ■ • C1F 


4.86 


4.21 


4.72 


4.21 


NH 3 ..-C1 2 


4.88 


4.(11 


5.02 


4.81 


H 2 0- • • C1F 


5.36 


5.21 


5.83 


5.18 


NH 3 - ■ • C1F 


10.62 


10.64 


11.07 


11.12 


MSE 




-0.24 


0.17 


-0.16 


MAE 




0.25 


0.33 


0.31 


Dipole-dipole interaction complexes 










H 2 S ■ ■ ■ H 2 S 


1.66 


1.44 


1.42 


1.54 


HC1- ■ • HC1 


2.01 


1.62 


1.71 


1.69 


H 2 S- ■ • HC1 


3.35 


3.26 


3.23 


3.38 


CH 3 C1- • ■ HC1 


3.55 


3.08 


3.24 


3.17 


HCN- ■ • CH 3 SH 


3.59 


3.68 


3.65 


3.72 


CH 3 SH- ■ • HC1 


4.16 


4.82 


4.84 


4.87 


MSE 




-0.07 


-0.04 


0.01 


MAE 




0.32 


0.28 


0.28 


Weak interaction complexes 










He- • ■ Nc 


0.04 


-0.02 


0.04 


0.00 


Ho- • • Ar 


0.06 


-0.03 


0.07 


0.00 


Nc- • ■ Nc 


0.08 


-0.03 


0.07 


-0.02 


Ne- ■ ■ Ar 


0.13 


-0.02 


0.13 


-0.01 


CH 4 - ■ • Nc 


0.22 


0.13 


0.20 


0.13 


C 6 H 6 ...Nc 


0.47 


0.44 


0.62 


0.24 


CH 4 -. CH 4 


0.51 


0.68 


0.51 


0.60 


MSE 




-0.05 


0.02 


-0.08 


MAE 




0.10 


0.03 


0.11 


Hydrogen-bonded DNA base pairs 










G- ■ ■ A HB 


-11.30 


-12.63 


-10.82 


-13.48 


C- • • G WC 


-30.70 


-31.42 


- -30.66 


-32.45 


G- • C WC 


-31.40 


-31.34 


-30.67 


-32.28 


MSE 




-0.67 


0.42 


-1.60 


MAE 




0.70 


0.42 


1.60 


Interstrand base pairs 










G- ■ • G IS 


-5.20 


-5.38 


-4.80 


-5.43 


G- ■ ■ G IS 


0.80 


0.88 


2.38 


1.15 


C • • C IS 


3.10 


3.34 


3.75 


3.29 


MSE 




0.04 


0.88 


0.10 


MAE 




0.17 


0.88 


0.25 


Stacked base pairs 










A- ■ • G S 


-6.50 


-6.84 


-3.61 


-6.66 


C • • G S 


-12.40 


-11.09 


-8.57 


-10.83 


G- ■ ■ C S 


-11.60 


-11.26 


-8.78 


-11.06 


MSE 




0.44 


3.18 


0.65 


MAE 




0.66 


3.18 


0.76 


MSE 




-0.11 


0.50 


-0.15 


MAE 




0.31 


0.61 


0.43 





